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Q ; ABSTRACT 

^ ■ We use a ray-tracing technique to compute the observed spectrum of a thin 

>- ■ accretion disk around a Kerr black hole. We include all relativistic effects such as 

frame-dragging, Doppler boost, gravitational redshift, and bending of light by the 

■ gravity of the black hole. We also include self-irradiation of the disk as a result of 

■ light deflection. Assuming that the disk emission is locally blackbody, we show 
how the observed spectrum depends on the spin of the black hole, the inclination 

■ of the disk, and the torque at the inner edge of the disk. We find that the effect of 
o a nonzero torque on the spectrum can, to a good approximation, be absorbed into 
q ■ a zero-torque model by adjusting the mass accretion rate and the normalization. 

We describe a computer model, called KERRBB, which we have developed for 
fitting the spectra of black hole X-ray binaries. Using KERRBB within the X- 
ray data reduction package XSPEC, and assuming a spectral hardening factor 
^ ■ /coi = 1.7, we analyze the spectra of three black hole X-ray binaries: 4U1543-47, 

XTE J1550-564, and GRO J1655-40. We estimate the spin parameters of the 
black holes in 4U1543-47 and GRO J1655-40 to be a/M ~ 0.6 and ~ 0.6 - 0.7, 
respectively. If / co i ~ 1.5 — 1.6, as in a recent study, then we find a/M ~ 
0.7 — 0.8 and ~ 0.8 — 0.9, respectively. These estimates are subject to additional 
uncertainties in the assumed black hole masses, distances and disk inclinations. 

Subject headings: black hole physics — accretion, accretion disks — radiation 
mechanisms: thermal — X-rays: binaries 
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1. Introduction 

Although the standard theory of thin accretion disks around black holes was developed 
over thirty years ago (Pringle & Rees 1972; Shakura & Sunyaev 1973; Novikov & Thorne 1973; 
Lynden-Bell & Pringle 1974), a straightforward confrontation of the model with observations 
is still not possible because of the challenging task of computing the observed spectrum of 
an accretion disk around a Kerr black hole. When fitting the soft X-ray spectra of black 
hole binaries, the multi-temperature disk model DISKBB is often used, which describes 
an approximate Newtonian model of a thin disk (Mitsuda et al. 1984; Makishima et al. 
1986). This model has been criticized by Gierlihski et al. (2001) for not applying the proper 
boundary condition at the inner edge of the disk. More seriously, the model does not include 
relativistic effects. 

Various efforts have been made to include relativistic effects and the appropriate bound- 
ary condition in calculating the spectra of black hole accretion disks (Cunningham 1975, 
1976; Hanawa 1989; Ebisawa et al. 1991; Zhang et al. 1997; Agol & Krolik 2000; Gierlihski 
et al. 2001; Ebisawa et al. 2003; Davis et al. 2004). However, accurate quantitative anal- 
ysis of observational data, e.g., trying to determine the radius of the disk inner edge from 
spectral data with a view to estimating the spin of the black hole, requires a sophisticated 
model that treats all relativistic effects including the self-irradiation of the disk caused by 
light deflection by the central black hole. At present, no such complete model exists that is 
suitable for inclusion with standard data reduction software, e.g., the X-ray spectral fitting 
package XSPEC (Arnaud 1996). As high-quality observational data on black-hole accretion 
disks become increasingly available (e.g., McClintock & Remillard 2004), the need for such a 
model has become urgent. In this paper we present the tools for a general relativistic model 
of an accretion disk around a Kerr black hole, and describe a numerical code for computing 
the observed blackbody spectrum of the disk. The code is efficient, and we have developed 
the appropriate software for interfacing it with XSPEC for analyzing spectral data. 

Two basic approaches have been described in the literature for calculating the observed 
spectrum of an accretion disk around a Kerr black hole. The first method makes use of 
"transfer functions" (Cunningham 1975, 1976; Laor 1991; Speith et al. 1995; Agol & Krolik 
2000; Dovciak et al. 2004). In this approach, all the information about the Doppler boost due 
to the disk rotation as well as the relativistic light deflection in the vicinity of the black hole 
is contained in a transfer function, which operates as an integration kernel for calculating 
the overall disk spectrum. The transfer function has been calculated and discussed in great 
detail by Speith et al. (1995), who have provided a fast and easy-to-use computer program 
to do the calculations. The second method uses the technique of "ray-tracing" (Rauch & 
Blandford 1994; Fanton et al. 1997; Cadez et al. 1998; Miiller & Camenzind 2004; Schnittman 
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& Bertschinger 2004). In this method one divides the image of the disk on the observer's 
sky into a number of small elements. For each image element, the orbit of a photon is traced 
backward from the observer by following the geodesies in a Kerr spacetime, until the orbit 
crosses the plane of the disk. The flux density of the radiation emitted by the disk at that 
point, as well as the redshift factor of the photon and the angle between the wavevector of 
the photon and the normal to the disk surface, are calculated. The observed flux density 
contributed by each image element is thus obtained, and summing over all the elements gives 
the total observed flux density of the disk (see eq. [E7] below). 

In this paper we use the ray-tracing approach, which we find more straightforward for 
numerical computations. We follow the procedures described by Fanton et al. (1997) and 
Cadez et al. (1998), who used elliptic integrals to simplify the calculation of the orbit of a 
photon in the background of a Kerr black hole (see also Rauch & Blandford 1994). Fanton et 
al.'s code was written for computing line emission from disks, whereas our code is designed 
for the continuum blackbody emission. We have made several improvements to Fanton et 
al.'s calculations; in particular, we include the effect of returning radiation, and we allow a 
nonzero torque to be set at the inner edge of the disk. 

Using our code, we have calculated spectra corresponding to a three-dimensional grid of 
models spanning different values of the black hole spin parameter a* = a/M, disk inclination 
angle ^obs, and dimensionless torque parameter 77 (defined in eq. [2]). We have also developed 
associated software which goes by the model name KERRBB for use with XSPEC. KERRBB 
reads in the above table of spectra and uses it to fit spectral data. In addition to the three 
parameters, a*, $ bs, and rj spanned by the table, other parameters such as the mass of the 
black hole M, the distance to the source D, the mass accretion rate of the disk M, and the 
spectral hardening factor / col (see eqs. [E5] and [E6]), are also included. The user can also 
decide whether to assume isotropic emission from the disk surface or to include a standard 
limb-darkening law (eq. [D20]). The model works efficiently within XSPEC and some sample 
results are presented. 

The paper is organized as follows. In §2 we summarize the assumptions behind our 
model, and in §3 we present the basic mathematical formalism. We discuss the effect of 
returning radiation in §3.1 and show examples of calculated spectra in §3.2, where we explain 
how the spectrum is affected by parameters such as the spin of the black hole, the inclination 
of the disk, and the torque on the inner edge of the disk. In §4 we apply KERRBB to spectral 
data on three black hole X-ray binaries, 4U1543-47, XTE J1550-564, and GRO J1655-40, 
and compare the results to those obtained with other models. We conclude with a summary 
in §5. 

The ray-tracing technique involves complicated mathematics to describe the orbits of 
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photons in Kerr spacetime. Although many of these formulae can be found in the literature 
(e.g., Chandrasekhar 1983; Rauch & Blandford 1994; Fanton et al. 1997; Cadez et al. 1998), 
we feel it is important to present the relevant mathematics in full detail to make the paper 
more useful to future workers in this field. To avoid distracting the reader, we give the 
technical details in Appendixes A-E. 

In Appendix F we compare KERRBB with GRAD — a subroutine in XSPEC for calcu- 
lating the blackbody spectrum of a Keplerian disk around a nonrotating (a* = 0) black hole. 
We show that KERRBB and GRAD give consistent results, once some errors in GRAD are 
corrected. 



2. Basic Assumptions 

Throughout the paper we use units in which G = c = h = 1, where G is the Newtonian 
constant, c is the speed of light, and h is the Planck constant. We use cylindrical coordinates 
(t,r,z,ip) as described by Page & Thorne (1974). In these coordinates, the rotation axis of 
the black hole is along the z-axis, and the equatorial plane corresponds to z — 0. 

We consider a geometrically thin and optically thick Keplerian accretion disk around a 
Kerr black hole, with the spin axis of the black hole perpendicular to the disk plane. The 
black hole has a mass M and a specific angular momentum a, where — M < a < M. The disk 
has an inner boundary of radius r in , and an outer boundary of radius r out 3> r in (KERRBB 
assumes r out = 10 6 M). 

In the standard theory of accretion disks, it is usually assumed that the torque at the 
inner boundary of the disk is zero. However, this assumption has only been justified for 
nonmagnetized or weakly magnetized flows (Muchotrzeb & Paczynski 1982; Abramowicz & 
Kato 1989) and for very thin disks (Afshordi & Paczynski 2003). Recent theoretical works 
on accretion disks have suggested that a nonzero torque at the inner boundary can arise from 
a magnetic field which either connects a disk to a central black hole, or couples a disk to the 
material in the plunging region (Krolik 1999; Gammie 1999; Agol & Krolik 2000; Li 2000, 
2002a,b, 2004; Wang et al. 2002, 2003; Uzdensky 2004a,b). Although this issue is still under 
debate (Paczynski 2000; Armitage et al. 2001; Hawley & Krolik 2002; Afshordi & Paczynski 
2003; Li 2003a,b), to make our model as general as possible we assume that the torque at 
the inner boundary of the disk can have any nonnegative value. 

In the presence of a torque g m > at the inner edge, the total power of the disk, i.e. 
the net amount of energy flowing out of the disk per unit time as measured by an observer 
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at infinity, is (Li 2002a) 1 

Aotal = #iiAn + C in M , (1) 

where fli n is the angular velocity of the disk inner boundary M is the mass accretion rate of 
the disk, and ei n = 1 — E- m is the specific gravitational binding energy at the inner boundary, 
where £j n is the specific energy of disk particles at the inner boundary. 

Equation (1) shows that the total power of the disk comes from two sources: the grav- 
itational binding energy between the disk and the black hole, and a contribution from the 
torque at the inner boundary of the disk (whenever the torque is nonzero). If the torque 
is produced by the black hole, then the power source for the second component is the spin 
energy of the black hole. 

Let us define a dimensionless parameter 

77 EE , (2) 

e in M 

which measures the ratio of the power from the torque to the power from the gravitational 
binding energy of the accreting gas. Then we have 

Aotal = (1 + VhnM . (3) 

In this paper we treat 77 as a free parameter that can have any nonnegative value. It is useful 
to define an effective mass accretion rate 

M cS = (l + V )M, (4) 

such that the total power of the disk is simply 

Aotal = CinMeff . (5) 



With our conventions, 77 = corresponds to the case when the torque at the inner 
boundary of the disk is zero (we call this the standard case), and 77 = 00 corresponds to 
M = (we call this the nonaccreting case). In these two extreme cases, all the power of the 
disk comes from either accretion (77 = 0) or from the torque at the inner boundary (77 = 00). 



^^Whcn the effect of returning radiation is considered, the total power of the disk corresponds to the 
energy-at-infinity carried away from the disk per unit coordinate time by the photons that permanently 
leave the disk — they either escape to infinity or fall into the black hole. See §3.1 for details. 
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The nonaccreting case is similar to the case of a "dead disk" around a magnetized 
spinning neutron star (Syunyaev & Shakura 1977), in the sense that both correspond to a 
zero mass accretion rate. However, the nonaccreting case considered in this paper is not 
necessarily in a "dead" state. If the black hole rotates fast and the inner boundary of a 
nonaccreting disk is located at the marginally stable orbit close to the horizon of the black 
hole, then a disk powered by the spin energy of the black hole can, in fact, be very bright 
(Li 2004). 

The model and the computer code described in this paper apply to a geometrically thin 
Keplerian disk with its inner edge located at any r; n > r ms , where r ms is the radius of the 
marginally stable orbit. For simplicity, in all plots presented in the paper, and in the model 
which we have developed for use with XSPEC, we assume that r in = r ms . 

We assume that the disk radiates like a blackbody. However, due to the complicated 
scattering processes in the disk atmosphere (predominantly electron scattering and Comp- 
tonization), the color temperature T co \ of the emitted radiation is generally higher than the 
effective temperature T eff of the disk (Ross et al. 1992; Shimura & Takahara 1993, 1995; 
Davis et al. 2004). To take this effect into account, we follow Shimura & Takahara (1995) 
and Ebisawa et al. (2003) and assume that the ratio / co i = T col /T eff is a constant (see eqs. 
[E4]-[E6]). We take / co i = 1.7, the mean value recommended by Shimura & Takahara (1995), 
when we model observations in §4. Undoubtedly, it is an over simplification to model the 
spectral modification due to Comptonization and electron scattering with a single scaling 
parameter, but this is a standard approach in the literature and is the best that we can do 
at this time. 



3. Mathematical Formalism for the Calculation of the Observed Spectrum of 

Disk Radiation 

3.1. Effects of the returning radiation 

Because of the gravity of the central black hole, not all of the radiation emitted by the 
disk escapes to infinity: a part of it, which we call "returning radiation," returns to the disk 
or is permanently captured by the black hole. The returning radiation that strikes the disk 
will interact with the disk particles and eventually be scattered or absorbed. For simplicity, 
we make the following assumption: All the radiation returning to the disk in the region 
beyond the inner boundary is absorbed by the disk, then reprocessed and reradiated. All 
the radiation returning in the region inside the inner boundary (i.e., the plunging region) is 
advected or scattered inward by the infalling gas (which has a large inward velocity), and 
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is captured by the black hole (Agol & Krolik 2000). In this subsection we study the effects 
of the returning radiation on the emission of the disk and the evolution of the central black 
hole. 

The mathematics for studying the effect of the returning radiation is presented in Ap- 
pendix D. In brief, considering the effect of returning radiation, at each point in the disk 
the net flux density F is composed of two nonnegative components: an outgoing component 
F out which represents the flux density of the energy emitted by the disk in situ; and an 
ingoing component F in , which is the flux density in the returning radiation, i.e. radiation 
emitted elsewhere and focused back onto the disk by the gravity of the black hole. In the 
steady state, the net flux density F is determined by the balance of the energy and angular 
momentum in the disk: F = F + F s , where F is the standard solution for the flux density 
when the effect of returning radiation is ignored (eq. [Dll]), and F s represents the work 
done by the returning radiation on the disk (eq. [D12]). The total outgoing flux density is 
then F out = F + F- m = F + F in + Fs (eq. [D17]). The inclusion of Fs in our energy balance 
equation is the main difference between our calculations and those of Agol & Krolik (2000). 
Note, however, that this has very little effect on the flux or spectrum of a standard disk. 

The self-irradiation of the disk arising from returning radiation is essentially a nonlocal 
process. Both F[ n and Fs of the incoming radiation — and also F out at the point on the disk 
where the incoming photons cross the disk — are functionals of F out at the point on the disk 
where the incoming photons were emitted. Therefore, F ont (and thus F; n and F$) must be 
obtained by solving the following functional equation 

-^out = Fq + -Fin [-Fout] + Fs[F out ] . 

Self-consistent solutions for F out , F in , and F s can be obtained by an iterative method using 
ray-tracing (see Appendix D). 

In Figure 1 we show the solutions of F in and Fs for the returning radiation as functions 
of disk radius, for the case of a Kerr black hole of a = 0.999M with rj = (upper panel) 
and r] = oo (lower panel). The radiation emitted by the disk is assumed to be isotropic in 
the disk frame, rather than limb-darkened (see eq. [D20] in Appendix D). As mentioned 
earlier, in all figures in this paper we assume that the inner boundary of the disk is at 
the marginally stable orbit, r- m = r ms . For comparison, the outgoing flux density when the 
returning radiation is ignored (Fq, eq. [Dll]) is also shown. The flux densities -F n and Fs 
always have the following asymptotic behaviors: F in oc r~ 3 , Fs oc r~ 7 / 2 for r ^> r in . Thus, 
asymptotically, F in behaves as if it is produced by a "lamp" on the axis of the black hole, 
while Fs behaves like it is produced by a torque at the inner boundary of the disk. We 
remind the reader of the asymptotic behavior of F : F oc r~ 3 when rj = 0, F oc r~ 7 / 2 
when i] = oo. This asymptotic behavior is the same as that in the Newtonian case (see, e.g., 
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Syunyaev & Shakura 1977), which is not surprising since at large radii the effect of relativity 
is unimportant. 

From Figure 1 we see that, for the case rj = (upper panel), the flux density of the 
returning radiation is always dominated by the original outgoing flux density F Q , except 
near the inner boundary of the disk where F and F$ approaches zero but F in remains finite. 
When t] is large (lower panel), at small radii the flux density is dominated by the original 
outgoing flux density F , while at large radii the flux density is dominated by that of the 
returning radiation (F m ). Since at large radii Fi n is always oc r~ 3 , which is the same as the 
asymptotic behavior of Fq(t] = 0), the lower panel of Figure 1 indicates that at large disk 
radii the signature of a large torque at the inner boundary is smeared out by the returning 
radiation. This makes it hard to detect the torque at the disk inner boundary by observing 
the spectrum of the disk (see §3.2 for more discussion). 

From Figure 1 we also see that F$ is always less important than F m . 

A photon emitted by the disk has three possible fates: it may be captured by the black 
hole, return to the disk, or escape to infinity. Thus, when the returning radiation is consid- 
ered, we have three different definitions for the "power" of the disk: 1) The total emission of 
the disk, which is the total energy- at- infinity emitted by the disk per unit coordinate time 
(i.e., the total energy-at-infinity carried away from the disk per unit coordinate time by all 
photons emitted by the disk, whether the photons escape to infinity, fall into the black hole, 
or return to the disk) 



where E' is the specific energy of a disk particle on a circular orbit of radius r (Page & 
Thorne 1974; Thorne 1974). 2) The total power of the disk, defined by equations (D22) and 
(1), which is the total energy-at-infinity carried away from the disk per unit coordinate time 
by the "net" radiation, i.e., the photons that permanently leave the disk — they either escape 
to infinity or fall into the black hole. 3) The total luminosity of the disk, which is the total 
energy-at-infinity carried away from the disk per unit coordinate time by the photons that 
escape to infinity. 

By definition, the total emission of the disk, P cm it, can be decomposed into three com- 
ponents: P cm it,BH! for the radiation that falls into the black hole; -P e mit,ret, f° r the radiation 
that returns to the disk; and P e mit,esc, for the radiation that escapes to infinity. Therefore, 
we can define three fractions 




(6) 



emit,BH 
p 

1 emit 



emit, ret 



emit, esc 



(7) 



^BH = 



esc — 



emit 



emit 
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Clearly, the three ratios satisfy lbu + <ret + W = 1- 



By the conservation of energy, we must have 



P 



emit, esc 



cmit,BH 



'total j 



(8) 



where £ to tai is the total power of the disk 



So we can define two other ratios 




cmit,BH 



esc — 



emit, esc 



(9) 



total 



■total 



which respectively represent the fraction of the energy going into the black hole and the 
fraction of the energy escaping to infinity in the "net" energy radiated by the disk. The two 
ratios must satisfy /bh + /esc = 1- 

The fractions with respect to the "total" radiation, t ret , £bh, and l csc , are shown as 
functions of the spin of the black hole in Figure 2. Again, we assume that the radiation 
emitted by the disk is isotropic in the disk frame. Two extreme cases are shown: a standard 
Keplerian disk, where g m = (equivalent to rj = 0, thin lines); and a nonaccreting disk, 
where M = but g m ^ (equivalent to rj = oo, thick lines). For the rj = case, we show 
the spin of the black hole from a = to a = 0.9999M. By considering the thermodynamics 
of black holes, Agol & Krolik (2000) argued that a nonaccreting disk (77 = 00) can exist 
only for a > 0.3594M. A physical explanation for this is that a black hole rotates faster 
than the inner boundary of the disk (at the marginally stable orbit) only if a > 0.3594M. 
Therefore, if a < 0.3594M a black hole cannot exert a positive torque on the disk (Li 2000, 
2002a). Hence, for the case of a nonaccreting disk, we show the spin of the black hole from 
a = 0.3594M to a = 0.9999M. 

The corresponding fractions with respect to the "net" radiation, / BH and / esc , are shown 
in Figure 3 (upper panel). To check the conservation of energy, i.e. equation (8), the 
difference between the computed /bh + /esc and 1 is also shown (lower panel). Within the 
errors of the computation the conservation of energy is confirmed. 

From Figures 2 and 3 we see that the effect of the returning radiation crucially depends 
on the spin of the black hole and the torque at the inner boundary of the disk. For a 
standard accretion disk (rj = 0) around a Schwarzschild black hole (a = 0), 1.7% of the total 
radiation emitted by the disk returns to the disk, 0.66% is captured by the black hole, and 
the remaining 97.6% escapes to infinity. For a standard accretion disk (77 = 0) around a Kerr 
black hole of a = 0.9999M, 27% of the total radiation emitted by the disk returns to the 
disk, 4% is captured by the black hole, and 69% escapes to infinity. For a nonaccreting disk 
(77 = 00) around a Kerr black hole of a = 0.9999M, 59% of the total radiation emitted by 
the disk returns to the disk, 7% is captured by the black hole, and 34% escapes to infinity. 
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Therefore, the effect of the returning radiation is most important for a fast spinning black 
hole and a directly rotating disk (rotating in the same direction as the black hole) with a 
large torque at its inner boundary. The obvious reason for this is that as a goes up the inner 
boundary of the disk shrinks, and when g in > more energy is dissipated in and radiated 
from the inner region of the disk, which is close to the central black hole. 

In terms of the "net" radiation of the disk, for a standard accretion disk around a 
Schwarzschild black hole, about 0.7% of the "net" radiation is captured by the black hole, 
and the remaining 99.3% escapes to infinity. For a standard accretion disk around a Kerr 
black hole of a = 0.9999M, about 6% of the "net" radiation is captured and 94% escapes 
to infinity. For a nonaccreting disk around a Kerr black hole of a = 0.9999M, about 17% of 
the "net" radiation is captured and 83% escapes to infinity. 

Similar results have been obtained by Agol & Krolik (2000) using the "transfer function" 
approach. However, in their calculations, Agol and Krolik ignored the stress of the returning 
radiation (i.e., the term F$) and obtained a fraction for the radiation captured by the black 
hole in the "net" radiation that is somewhat smaller than the fraction that we have obtained. 
For example, for a nonaccreting disk around a Kerr black hole of a = 0.9999M, they find 
/bh = 15%, while we obtain / B h = 17%. 

Figures 1-3 show that when the black hole is rotating rapidly and there is a large torque 
at the inner boundary of the disk the effect of returning radiation is extremely important. 

In order to study the effect of the radiation that is captured by the black hole on the 
spinup/spindown of the black hole, following Thorne (1974) we define a capture function C 
for each photon emitted by the disk at radius r: C — 1 if the photon is eventually captured 
by the black hole, and C — if the photon escapes to infinity or returns to the disk. The 
calculation of the capture function is described in Appendix A. 3. 

Here we focus on the case when the stress at the inner boundary of the disk is zero, 
i.e., <7i n = 0. Then, the torque that spins up/down the black hole comes from two sources: 
gas accreted from the disk, and radiation captured by the black hole. Then, the resultant 
spinup/spindown of the black hole is governed by (Thorne 1974) 

-p-= ft + C * -2a., o.^, (10) 
dlnM M (Et + C E ) 

where L\ is the specific angular momentum of disk particles at the inner boundary, and 



Ce = 



2vrr2 



CT(-n t ) cos 9dQ 



n 4 



foutrdr , (11) 
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Cl = ^/"(^ CTnyCoaOdn) f^rdr , (12) 

where r g = M is the gravitational radius of the black hole, n t and are given by equa- 
tion (D5). See Appendix D for the meanings of 9, fi, Q + , and For convenience, in 
equations (11) and (12) we have used a dimensionless outgoing flux function / out defined by 

-fout = — 5" /out • (13) 

The function da*/d\nM, which is defined by equation (10) and is a function of a* only 
when r in = r ms , is plotted in Figure 4 (upper panel). The value of a at which da*/dlnM = 
gives the equilibrium spin of the black hole [which is called the "canonical state" by Thorne 
(1974)]: a = a eq . When the black hole is in a state of a < a eq , the effect of accretion from 
the disk dominates, which will spin up the black hole, until the equilibrium state a = a cq 
is reached. On the other hand, when the black hole is in a state of a > a eq , the effect of 
the capture of photons dominates, which will spin down the black hole to the equilibrium 
state a = a eq . When the disk emission is isotropic, we obtain a eq = 0.9983M, and when 
the disk emission is limb-darkened, we obtain a eq = 0.9986M. 2 These results show that the 
limb-darkening effect does not significantly affect the gross disk radiation. 

The efficiency of an accretion disk in converting rest mass into outgoing radiation is 
defined by the ratio of the total luminosity of the disk to the mass-energy accretion rate as 
measured at infinity, which is given by 1 — Ej n — (e (Thorne 1974). In the lower panel of 
Figure 4 we show the efficiency of a standard accretion disk as a function of the spin of the 
black hole when the effect of the returning radiation is considered. In the canonical state, 
the total efficiency is 0.309 when the disk emission is isotropic and 0.315 when the emission is 
limb-darkened. When the effect of the returning radiation is ignored, the radiation efficiency 
of a standard disk is 0.326 when a = 0.9983M, and 0.331 when a = 0.9986M. 

For comparison, the results of Thorne (1974) are also shown in the lower panel of Figure 4 
(the two plus signs). Without considering the effect of the radiation returning to the disk, 
Thorne obtained: a cq = 0.9978M when the disk radiation is isotropic, and a cq = 0.9982 
when the disk radiation is limb-darkened. The corresponding radiation efficiencies of the 
disk are respectively 0.302 and 0.308. 



2 These numerical results are correct only if accretion from a thin disk and capture of photons is the 
unique process to spin up/down the black hole, and if the torque at the inner boundary of the disk is zero. 
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3.2. Blackbody radiation spectrum from a Keplerian accretion disk 

The energy flux density of the blackbody radiation emitted by an accretion disk around 
a black hole as observed by a remote observer — FE obs , where E ^, s is the photon energy — is 
given by equation (E7). However, for the purpose of comparison with X-ray observations, 
it is more convenient to use the photon number flux density rather than the energy flux 
density. So we define N ohs = F Eobs /E ohs . 

For convenience, let us define 

d&obs = (JpJ dD, ohs , (14) 

where D is the distance from the observer to the black hole and dil b s is the element of 
the solid angle subtended by the image of the disk on the observer's sky. With the above 
definition, <if2 bs defined above is independent of the distance D. Then, by equations (E7) 
and (E4), the formula for calculating N ohs is 



Eobs\ f G^obs 



cxp 



,. (Mobs.) a-if- 1 / 4 
P I keV ) y J out 



- 1 



where g is the photon redshift (defined by eq. [E3]), / ou t is defined by equation (13), and 

= 0.07205 /-f (J^ (k^) Photons keV^^m^ 2 sec" 1 , (16) 
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where / co i is the spectrum hardening factor defined by equation (E5). 
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Some examples of the observed spectra calculated with our ray-tracing code are shown in 
Figures 5-9. In Figures 5-8 we assume the disk radiation is isotropic in the frame corotating 
with the disk, i.e., there is no limb-darkening, whereas in Figure 9 we include limb-darkening 
according to equation (D20). 

In the upper panel of Figure 5, we show the dependence of the observed spectrum 
on the spin of the black hole. The disk has i] = (zero-torque) and an inclination angle 
^obs = 30° (for other parameters see the caption of the figure). Different lines correspond to 
different spins of the black hole: a/M = 0, 0.5, 0.9, and 0.999 (left to right). We see that the 
spectrum becomes harder as the spin of the black hole goes up. Physically, this is caused 
by the fact that when M is fixed and a/M increases, the radius of the inner edge of the 
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disk (by assumption located at the marginally stable orbit) decreases, so that the disk has 
a higher radiation efficiency (see Fig. 4, lower panel) and a higher temperature. We also 
see another effect: as, a/M increases, the flux density at low energies increases even though 
this radiation comes from large radii where the spin of the black hole should have negligible 
effect. This increase in flux density is caused by the returning radiation. As a/M goes up, 
more radiation emitted by the disk in the inner region is focused back to the disk by the 
gravity of the black hole (Figs. 2 and 3). This increases the radiation of the disk even at 
large radii (corresponding to the low energy end of the spectrum). 

In the lower panel of Figure 5, we test the effect on the observed spectrum of the 
inclination angle of the disk. The case considered is for rj = and a/M = 0.9. Different 
lines correspond to different inclination angles of the disk: -# obs = 0°,40°,70° and 85°. At 
the low energy end, the flux density goes down as -i9 ob s increases. This is caused by the 
projection effect. The low energy radiation is primarily emitted by the disk at large radii, 
where the effect of relativity is not important. The projection causes the flux density of the 
disk radiation to be proportional to cos^obs- At the high energy end, the flux density goes 
up as $obs increases. This is caused by the effects of Doppler beaming and gravitational 
focusing. The high energy radiation is primarily emitted by the disk in the region near 
the black hole, where the orbital velocity of the disk is mildly relativistic so that special 
relativistic beaming boosts the disk radiation to higher energy. In addition, near the black 
hole, the gravity of the black hole is strong and focuses the disk radiation back to the disk 
plane, thereby modifying the projection effect. The joint action of the two effects leads to 
an enhancement at the high energy end of the observed spectrum. 

Figure 6 shows the effect of the returning radiation on the observed spectrum of the 
disk. The three panels correspond to rj = (upper panel, the standard disk case), rj — 1 
(middle panel), and rj = oo (lower panel, the nonaccreting case). The black hole has a spin 
a/M = 0.999, and the disk has an inclination angle -# obs = 30° (other parameters are given in 
the caption of the figure). The spectra when returning radiation is included are shown with 
solid lines, and those without the returning radiation are shown with dashed lines. We see 
that the returning radiation enhances the disk radiation (especially at the high energy end), 
and the effect is more prominent for a disk with a larger torque at its inner boundary. For the 
standard disk with 77 = (upper panel), where the power of the disk comes purely from disk 
accretion, the effect of the returning radiation is almost indistinguishable from the effect of 
a change in the mass accretion rate: The dotted line (almost coincident with the solid line) 
represents a disk spectrum minus the returning radiation with the same parameters except 
that the mass accretion rate is larger by a factor of 1.23. For the case rj = 1 (middle panel), 
for which the power of the disk comes equally from disk accretion and a torque at the inner 
boundary of the disk, the effect of the returning radiation can again be well approximated by 



-14- 



adjusting the effective mass accretion rate. The dotted line represents the spectrum without 
the returning radiation of a disk with the same parameters except that the effective mass 
accretion rate is larger by a factor of 1.7. In this case also, the dotted line agrees very 
well with the solid line. However, when 77 is very large, e.g. for the nonaccreting case with 
77 = 00 (lower panel), where the power of the disk comes purely from the torque at the inner 
boundary, the effect of the returning radiation is so prominent that it cannot be fitted by 
simply modifying the effective mass accretion rate. 

Figure 7 shows the effect of a nonzero torque at the inner boundary of the disk for the 
case of a Schwarzschild black hole (a = 0). Each panel corresponds to a different value of 
77: 0.1, 0.3, 1 and 16.49. The last value corresponds to the case that the total efficiency of 
the disk is equal to unity — the upper limit for the efficiency of a Keplerian disk around a 
nonrotating black hole. In each panel, the solid line is the spectrum when the disk has a 
nonzero torque at the inner boundary, and the dashed line is the spectrum when the disk has 
a zero torque with other parameters remaining the same. In particular, for both the solid 
and the dashed lines the effective mass accretion rate is M e g = 10 19 gsec _1 , which means 
that the two disks have the same total power (but different mass accretion rates M, see 
eqs. [4] and [5]). From the figure we see that the effect of the torque at the inner boundary 
is to make the spectrum harder. When the torque at the inner boundary is positive, more 
energy is dissipated and more radiation is emitted in the inner region of the disk where 
the temperature is higher. Also, the disk rotates faster there, so the Lorentz boost is more 
prominent. Both effects cause a hardening of the spectrum. 

In the case of the returning radiation, we showed earlier that its effect can be modeled 
very well by adjusting the mass accretion rate in a model without returning radiation. Can 
the effect of a nonzero torque be similarly absorbed by adjusting the parameters of a zero- 
torque model? The first panel of Figure 7 shows that for 77 = 0.1, which corresponds to 
a weak torque, a zero-torque model with the same power (M eS ) as the finite-torque model 
gives an almost indistinguishable spectrum. So, in this case, the answer to our question is a 
definite yes. For the cases of 77 = 0.3, 1, and 16.49, we see that just keeping M e g the same is 
not enough, since the dashed lines are noticeably different from the solid lines. However, by 
adjusting both M cff and / col of a zero-torque model, we can get a very good fit to the finite- 
torque models with 77 = 0.3 and 1, as shown by the dotted lines; the corresponding values of 
these parameters are given in the caption to the Figure. Only in the extreme case 77 = 16.49 
are we unable to fit the spectrum by adjusting M e g and / co i. Note that adjusting / co i is 
equivalent to adjusting the normalization, which could be done equally well by adjusting M, 
D, or cos79 bs- Since none of these parameters is known precisely in a real system, one always 
has some freedom in the normalization. These results indicate that for a Schwarzschild black 
hole the spectrum of a disk with a modest torque at the inner boundary can be fitted with 
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a zero torque at the inner boundary by adjusting M and the normalization. 

Similar results for the case of a Kerr black hole with a/M = 0.999 are shown in Figure 8. 
Now we find that the spectrum can always be fitted by a disk with a zero torque at the inner 
boundary by adjusting M e ff and / co i; in fact, we can do this even when r\ = oo. Interestingly 
the effect of a nonzero torque is more important for a non-rotating black hole than for 
a rapidly rotating hole, e.g., note that the spectrum is significantly hardened in Figure 7 
whereas there is almost no effect in Figure 8. This is caused by the effect of the returning 
radiation. For a nonrotating black hole, the effect of the returning radiation is not important 
due to the fact that the inner edge of the disk has a large radius, so the torque at the inner 
boundary produces a radiation flux density (~ r~ 7 / 2 at large radii) that is distinctly different 
from the radiation flux density arising from accretion (~ r~ 3 at large radii). For a fast 
rotating black hole, the effect of the returning radiation is important, it makes the radiation 
flux density of the disk more or less similar to that arising from accretion (both going as 
~ r~ 3 at large radii) (see Figs. 1-2 and the relevant discussions in §3.1). 

Figure 8 shows another effect of the torque at the inner boundary for a rapidly spinning 
hole: when the torque changes from zero to nonzero and other parameters (including M e ff) 
are left unchanged, the flux at infinity decreases (compare the solid and dashed lines). This 
is caused by the fact that when the torque at the inner boundary is nonzero more energy is 
dissipated in the inner region of the disk; for a rapidly rotating black hole the inner boundary 
of the disk is closer to the horizon of the black hole and so more radiation is focused to the 
equatorial plane. 

Overall, the results shown in Figures 6-8 suggest that when modeling observational 
data the effects of returning radiation and nonzero torque can be ignored since they can 
be absorbed by modifying the mass accretion rate and the spectral hardening factor of a 
zero-torque model without returning radiation. 

Finally, in Figure 9 we show the effect of limb-darkening (eq. [D20], Appendix D) on the 
observed spectra of the disk. As can be expected from equation (D20), compared to the case 
when the disk emission is isotropic, when the disk emission is limb-darkened we see more 
radiation when the disk has a low inclination angle, and we see less radiation when the disk 
has a high inclination angle. The effect of limb-darkening is most important when the disk 
is nearly edge-on: in this case the effect cannot be absorbed by simply adjusting the mass 
accretion rate and the spectral hardening factor. 
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4. Modeling the Spectra of Black Hole X-ray Binaries 

To test the performance of our model KERRBB and to compare it with other models of 
disk spectra, we have analyzed spectral data on three black hole X-ray binaries: 4U1543-47 
(hereafter U1543), XTE J1550-564 (hereafter J1550), and GRO J1655-40 (hereafter J1655). 
The data were obtained using the Proportional Counter Array aboard the Rossi X-ray Tim- 
ing Explorer (RXTE). The data themselves and the methods of data analysis are thoroughly 
described in Zimmerman et al. 2004 (hereafter Z04). These three black hole binaries have all 
had recent outbursts that have been analyzed in detail using RXTE observations. Park et 
al. (2004) studied the 2002 outburst of U1543; Sobczak et al. (2000) analyzed the outburst of 
J1550 in 1998-1999; and Sobczak et al. (1999) analyzed J1655's 1996-1997 outburst. In fit- 
ting the spectra of these sources, we followed as closely as possible the procedures described 
in these papers. 

In addition to KERRBB and the other disk models discussed below, we used four other 
supplementary XSPEC models in our fits. Primary among these was the power-law model, 
which has two parameters: the photon index, T, and a normalization constant, which we will 
call -KpL- The remaining three spectral components, which model the effects of interstellar 
absorption, a smeared Fe absorption edge, and a Gaussian Fe line, are relatively unimportant. 
Because our focus was on the disk models, we froze a number of the extraneous parameters 
in these other spectral components to minimize their influence on our fits by fixing the 
parameters at average values obtained from the previously cited papers (see Z04). The 
spectra were fitted using XSPEC over the energy range ^3 — 20 keV. All of the data sets 
correspond to the high/soft state, for which the accretion flow is believed to be dominated 
by a geometrically thin, optically thick disk component of emission that contributes > 90% 
of the flux (McClintock & Remillard 2004). For further details on the data and the analysis 
techniques, see Z04. 

In addition to (i) KERRBB, the other XSPEC models we considered were: (ii) GRAD: 
a general relativistic code for non-spinning black holes (Hanawa 1989); (iii) DISKPN: a disk 
model employing a pseudo- Newtonian potential (Gierlihski et al. 1999); (iv) EZDISKBB: 
a Newtonian model with a zero-torque inner boundary condition (Z04); and (v) DISKBB: 
a Newtonian model (Mitsuda et al. 1984) which includes a finite torque at the inner edge 
of the disk (Gierlihski et al. 1999; Z04). Note that KERRBB is the most complete model 
of these five since it includes the effects of general relativity and can handle a black hole 
with any spin and any nonnegative torque at the inner edge of the disk. In this section 
we assume that the torque is zero (rj = 0). Likewise the torque at the inner disk edge 
is zero for all of the other models except DISKBB (see below). For the specific case of a 
non-spinning black hole (a = 0), KERRBB should agree with GRAD, but there are in fact 
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some deviations because of an error in GRAD discussed in Appendix F. In order to facilitate 
comparison between KERRBB and GRAD, we did not include the returning radiation in 
KERRBB in the calculations reported here. As we have shown in §3.2, the effect is in 
any event small for 77 = 0, and can be absorbed in the fitted value of M. Of the other 
three models, DISKPN should be closest to KERRBB since it attempts to include some 
relativistic effects through the pseudo-Newtonian potential of Paczyhski & Wiita (1980). 
The main difference between the final two Newtonian models is that EZDISKBB includes a 
zero-torque boundary condition at the inner edge of the disk whereas DISKBB has a finite 
torque. As argued in Z04, the zero-torque condition is expected to be valid in a number of 
situations. Furthermore, even in those cases in which a finite torque might be expected, it 
is unlikely that the specific magnitude of the torque would be equal to the value assumed in 
DISKBB. 



4.1. Comparison of normalizations 

For our first comparison, we analyzed data from ten epochs each on U1543 and J1550 
(see Z04 for details of the particular observations). We fixed the mass of the black hole and 
the inclination of the disk at their respective estimated values of 9.4M and 20.7° for U1543 
(Orosz et al. 1998; Orosz 2004) and 1O.6M and 73.5° for J1550 (Orosz et al. 2002). We 
also assumed that the two black holes are not spinning (a = 0) and fixed the inner edge 
of the disk at the marginally stable orbit (r- m = r ms = 6M). The constancy of this inner 
disk radius over a wide range of X-ray luminosity has been established in synoptic studies of 
several black hole binaries (Tanaka & Lewin 1995; Sobczak et al. 2000). Although there are 
reasonable estimates of the distances to the sources, we left the distance D as an adjustable 
parameter; in effect, this parameter played the role of a normalization in these calculations. 
We also left the mass accretion rate M as a free parameter, and used a spectral hardening 
factor /coi = 1.7 (Shimura & Takahara 1995). 

Figure 10 shows the results of fitting the data on U1543 with the five models. We see 
a remarkably consistent pattern. KERRBB and GRAD agree almost perfectly in their esti- 
mates of both D and M; DISKPN has modest deviations from these two models; EZDISKBB 
deviates somewhat more; and DISKBB deviates in both parameters by an enormous factor. 
Table 1 gives the average ratio by which the estimated parameters obtained with each of 
the models deviates from the value obtained with KERRBB. Here and in what follows, we 
assume that the result from KERRBB is correct and view any deviation from it as a measure 
of the error in a particular model. 

To understand the patterns seen in Figure 10, we note that each model adjusts two 
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parameters to the data. These parameters are obtained essentially by fitting the integrated 
flux from the source and the position of the peak in the spectrum; the latter in effect measures 
the peak temperature (as observed at infinity) of the emission from the disk. For a non- 
rotating black hole, the radiative efficiency of KERRBB and GRAD is e = 1 — (8/9) 1 / 2 = 
0.057; that is, for a given M, the luminosity is equal to 0.057Mc 2 . In contrast, DISKPN has 
e = 0.0625 and EZDISKBB has e = 0.0833. Thus, these models produce the same luminosity 
with a smaller M. Moreover, for a given luminosity, the fully relativistic models KERRBB 
and GRAD include Doppler and gravitational redshift factors, whereas the other models do 
not. These cause additional deviations. The net effect is that the estimate of M obtained 
with DISKPN is lower than the correct relativistic result by a factor of ~ 0.6, while M with 
EZDISKBB is only a quarter of the correct value. The variations in the derived values of D 
with the various models are straightforward to understand — each model adjusts D so as to 
fit the observed flux oc luminosity//} 2 . 

Compared to the other models, DISKBB makes an extraordinarily large error: the 
estimate of M is reduced by more than a factor of 50 relative to KERRBB. The reason for 
this is well understood (e.g., Kubota et al. 1998; Gierlihski et al. 2001; Z04), though not 
widely appreciated. DISKBB assumes a finite torque at the inner edge and so its radiative 
efficiency is three times that of EZDISKBB: e = 0.25. Furthermore, this model has its 
temperature maximum at the inner edge of the disk, whereas all the other models have their 
maxima at larger radii. Therefore, it predicts a substantially larger value of T max for a given 
luminosity. Since the spectral fit tries to reproduce the position of the peak in the observed 
spectrum, both effects act in the same direction and lead to a large decrease in the estimated 
value of M. Despite this well-known deficiency in DISKBB, the model is still widely used 
to model spectra of black hole X-ray binaries. 

Some authors have used DISKBB but included correction factors to make the results 
more consistent with other models. For example, Zhang et al. (1997) and Gierlihski & Done 
(2004) apply correction factors to the temperature and the flux derived from DISKBB. Other 
authors have attempted to include these corrections in the spectral hardening factor / co i (see 
Davis et al. 2004 for a review of this topic). While these approaches have some merit, we feel 
it is better to use the correct model, viz., EZDISKBB if one wishes to assume a Newtonian 
model, DISKPN if one is interested in a pseudo-Newtonian potential, and KERRBB if one 
would like to include all relativistic effects including returning radiation, and would like to 
consider different black hole spins and inner torques. (GRAD is valid only for a non-rotating 
hole.) Indeed, if we apply the approach of Zhang et al. (1997) and Gierlihski & Done (2004) 
to DISKBB or EZDISKBB for U1543, where we expect M/D 2 oc l/g GR where g GR is the 
relativistic correction factor introduced by Zhang et al. (1997), we get M/D 2 ps 0.22 for the 
DISKBB results, and M/D 2 ^ 0.68 for the EZDISKBB results in Table 1. These corrections 
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improve somewhat the results of the models without relativistic corrections (M/D 2 0.18 
for DISKBB and 0.56 for EZDISKBB), but are still far from the correct relativistic results 
of KERRBB, which by definition imply M/D 2 = 1. 

In the case of J1550, which is nearly edge-on, we find that the relativistic corrections 
introduced by Zhang et al. (1997) do not improve the results of DISKBB and EZDISKBB 
at all. Indeed, from Table 1 of Zhang et al. (1997), by interpolation we have gcR ~ 114 for 
6 = 73.5° and a* = 0, which brings the results of DISKBB and EZDISKBB farther from the 
results of KERRBB (see Table 2). 

We turn now to a comparison of the results for J1550 obtained using the five models; 
these results are summarized in Figure 11 and Table 2. Similar patterns are seen as in 
the case of U1543, but there are also differences. The latter are all caused by the fact 
that U1543 is a nearly face-on system ("$ bs = 20.7°) whereas J1550 is a nearly edge-on 
system (^obs = 73.5°). For the three nonrelativistic models, viz., DISKPN, EZDISKBB, and 
DISKBB, the effect of inclination is straightforward: the flux simply decreases by a factor of 
cos$ b s , and there is no change in the spectral shape. However, the two relativistic models 
KERRBB and GRAD have additional effects. First, they have stronger Doppler beaming 
for an edge-on system like J1550 and consequently appear brighter and hotter for a given 
M. Second, because of light deflection, the effective projected area of the disk as viewed 
by the observer is larger than one might expect with a simple cos-$ bs scaling, especially 
for the hot inner regions of the disk. For both reasons, these two models are able to fit a 
given observed flux with a smaller value of M than one might expect by naively scaling from 
a face-on system like U1543. As a result the ratios in Table 2 corresponding to DISKPN, 
EZDISKBB, and DISKBB are higher by a factor of ~ 2 - 3 relative to those in Table 1. In 
addition, we see that the values obtained with GRAD are higher than those from KERRBB. 
This is the result of the error in GRAD pointed out in Appendix F. The error has almost 
no effect for a face-on system like U1543, but becomes more serious for edge-on systems. 

Finally, we note that the x 2 values of the fits vary somewhat erratically from one epoch 
to the next. Epochs 2 and 3 of U1543 are relatively bad for all models, while epochs 6 and 9 
of J1550 are bad for three of the models (DISKPN, EZDISKBB, and DISKBB) but perfectly 
reasonable for the other two models (KERRBB and GRAD). In the case of our third source, 
J1655, we obtained large x 2 values for all except the first two epochs. For this reason, we 
did not include this source for the calculations discussed in this subsection. 
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4.2. Estimating the black hole spin 

In the previous subsection, we treated the distance D to each source as a free parameter 
(equivalent to the normalization) and fitted D and M from the data. However, we do have 
independent estimates of the distances: 7.5 ± 1.0 kpc for U1543 (Orosz et al. 1998, 2002); 
5.91^ kpc for J1550 (Orosz et al. 2002); and 3.2 ± 0.2 kpc for J1655 (Hjellming & Rupen 
1995). By using this additional information we should in principle be able to constrain 
the spin of the black hole (Zhang et al. 1997). (Note that we fixed a = in §4.1.) In 
the case of KERRBB, the calculation is straightforward — we simply fix the value of D and 
let a* = a/M and M be the free parameters. For the other models, however, a* is not an 
adjustable parameter. In the case of the three non-relativistic models, DISKPN, EZDISKBB, 
and DISKBB, we allow the radius of the inner edge r- m and M to be the free parameters. 
Having fitted r in from the data, we then use the ratio ri n /M to estimate a*, assuming that the 
inner edge is at the marginally stable orbit corresponding to the particular spin parameter. 
This procedure is somewhat arbitrary, but is fairly standard practice in the literature. We 
have not attempted a similar exercise with GRAD. 

Figure 12 shows the estimated spin parameter a* for each of the 10 observations of U1543 
and J1550, and the first two observations of J1655 (with M = 7.OM and $ bs = 69.5°; Orosz 
& Bailyn 1997). We have ignored the remaining 8 observations of J1655 because (for some 
unexplained reason) the x 2 values are very large. In the case of U1543, KERRBB gives a 
consistent estimate of a* ~ 0.6 for all the epochs. The near-constancy of the estimate is 
notable, especially since the luminosity and the mass accretion rate do vary from one epoch 
to the next (see Fig. 10). The models DISKPN and EZDISKBB give different values of a*, 
suggesting that it is dangerous to use these models to estimate the black hole spin. The 
model DISKBB is especially poor — in the case of both U1543 and J1550, it gives estimates 
of the disk inner edge that are too large (> 9M) to be consistent with any choice of the 
spin parameter. The results are more variable in the case of J1550. We find with KERRBB 
values of a* ranging all the way from —0.8 to —0.1. In part this is because the spectra are 
less sensitive to the value of a* when the parameter is negative, but in part it might also 
reflect the quality of the data (e.g., the very uncertain distance to the source; see below). In 
the case of J1655, the two epochs that we have analyzed give a* ~ 0.6 — 0.7, slightly larger 
than for U1543. 

The ability to estimate a* depends on having independent estimates of the disk inclina- 
tion $obs, the black hole mass M, the distance to the source D, and the spectral hardening 
factor /coi. However, even for the best sources, there are substantial uncertainties in these 
parameters. Table 3 shows for one epoch of U1543 how the fitted value of a* varies as each 
of the four input parameters is allowed to range over its la range of uncertainty. Leaving 
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aside / co i, we see that the uncertainties in the other three parameters are not particularly 
severe for this favorable case of U1543. Even allowing for the uncertainties it appears that 
one could infer that a* > 0.5. The spectral hardening factor, however, introduces a large 
uncertainty. The value of this parameter is not well constrained. Shimura & Takahara (1995) 
suggested a value of 1.7, but said that the value could be anywhere in the range from 1.5 to 
1.9 (the range covered in Table 3). A recent comprehensive analysis by Davis et al. (2004), 
which includes better opacities, comes down in favor of somewhat smaller values in the range 
/coi ~ 1.5 — 1.6, depending on the luminosity of the disk. For such values, the spin estimates 
go up by a modest amount, to a* ~ 0.7 — 0.8 for U1543 and ~ 0.8 for J1655. Apart from 
the question of what is the correct value of f co \ to use, a more serious concern is that it is a 
rather severe simplification to model the spectral modification due to Comptonization and 
opacity effects with a single scaling parameter / co i. It is hard to quantify the error from this 
approximation. In our opinion, if one wishes to estimate black hole spin via spectral fitting, 
one must first develop more reliable spectral models of disk atmospheres. The work of Davis 
et al. (2004) is a first step in this direction. 

In the case of J1550, the distance is highly uncertain: D = 2.8 — 7.6 kpc (Orosz et al. 
2002). Correspondingly, a* is poorly constrained. We analyzed the 10 observations of this 
source using the two extreme la values of D and obtained values of a* ranging all the way 
from —1 to +0.7. This emphasizes once again that we cannot hope to constrain the spin of 
a black hole via spectral fitting unless we have accurate input parameters. 

5. Summary 

We have developed a ray-tracing computer code to calculate the spectrum of a thin 
accretion disk around a black hole, assuming that the emission from each point on the 
surface of the disk is locally blackbody-like with a constant spectral hardening factor / co i. 
The code includes all relativistic effects. It also includes the effect of self-irradiation whereby 
radiation emitted at one point on the disk is deflected by the gravity of the black hole and 
illuminates another part of the disk. The code can handle any value of the black hole spin, 
the disk inclination, the disk inner radius (r in > r ms ), and the torque on the inner edge (the 
dimensionless parameter 77). It also allows the user to choose between isotropic emission and 
a Chandrasekhar (1950) limb-darkening law (eq. [D20]). 

In terms of the scope of the calculations, the present work is not markedly different 
from previous studies. For instance, Fabian et al. (1989) and Stella (1990) calculated line 
emission from a relativistic disk around a non-rotating black hole, while Laor (1991) repeated 
the calculation for a rapidly rotating black hole with a specific value of a/M = 0.998 (i.e., a 
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nearly maximally rotating hole). Our code calculates the continuum emission, rather than 
the line emission, and is more comparable to the blackbody continuum models calculated 
by Hanawa (1989) and Ebisawa et al. (1991) for a disk around a non-rotating black hole 
and by Cunningham (1975, 1976) and Gierlihski et al. (2001) for a disk around a rotating 
black hole. Our model supersedes Cunningham's model by including a torque at the inner 
boundary of the disk, and that of Gierlihski et al. by including both the torque and the 
effect of the returning radiation. The calculations of Agol & Krolik (2000) are very close 
to the present work. The only difference is that they did not include the additional energy 
released as a result of the torque applied by the returning radiation (the term Fs in §3.1 and 
Appendix D). As we have shown, this term is quantitatively not important. 

Using the ray-tracing code, we have calculated tables of spectra and developed a model 
called KERRBB for use with the X-ray spectral fitting package XSPEC. This model can 
handle black hole spins in the range a* = a/M = —1 to +0.9999, inclinations in the range 
"#obs = 0° to 85°, and torque parameters in the range rj — to 1. One can choose either 
isotropic emission or limb-darkened emission from the disk (eq. [D20]). In addition, as usual, 
one can set (or fit) the black hole mass M, the distance to the source D, the mass accretion 
rate M, and the spectral hardening factor / co i. KERRBB thus goes well beyond other 
models that are presently available in XSPEC for fitting the continuum emission from thin 
disks. The other models are based on a Newtonian potential (DISKBB, EZDISKBB) or a 
pseudo- Newtonian potential (DISKPN), or if they solve the relativistic problem it is only 
for a specific case, e.g., a non-rotating black hole (GRAD, see Appendix F). KERRBB is 
more general and complete than any of these models. One limitation of KERRBB is that the 
present implementation in XSPEC only considers disks that extend down to the marginally 
stable orbit (r- m = r ms ). Disks that evaporate into a corona/ADAF at radii r- m > r ms are not 
yet modeled, but could easily be incorporated in the future if there is a need. 

The effects of some of the model parameters on the observed spectrum are illustrated in 
Figures 5-9. The black hole spin parameter a/M and the inclination angle -$ bs are seen to 
have a strong effect on the spectrum. These model parameters must clearly be included when 
fitting spectral data. The effect of the returning radiation appears to be less important. At 
least for torque parameters rj < 1, its effect can be almost completely absorbed by a rescaling 
of the mass accretion rate. Likewise, the effect of limb-darkening, which is important only 
when the disk is nearly edge-on, can otherwise be absorbed by adjusting the mass accretion 
rate and the normalization (Fig. 9). 

Afshordi & Paczyhski (2003) have argued that a thin accretion disk will have a negligible 
torque at its inner edge. Magnetohydrodynamic simulations of disks have found nonzero 
torques (Hawley & Krolik 2002), but the magnitude appears to be relatively small (77 ~ 0.2) 
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even for relatively thick disks (H/R ~ 0.1 — 0.2). For such values of 77, and indeed even 
for larger values up to 77 ~ 1, Figures 7 and 8 show that the effect of a finite torque can 
be modeled with a zero-torque disk by simply adjusting the mass accretion rate and the 
normalization (however, see the caveat in the following paragraph). Since M is almost 
always fitted to the data, and the normalization is rarely known accurately a priori because 
of uncertainties in the black hole mass, source distance, disk inclination, etc., it appears that 
there is no particular advantage to including a finite torque in fitting spectral observations. 
Nevertheless, we have retained rj as a parameter in KERRBB for the occasions when it might 
be needed. 

The following additional caveat must be noted with respect to models with a finite 
torque at the disk inner edge. Most such models have continued viscous dissipation at radii 
inside r ms (e.g., Krolik 1999; Gammie 1999). KERRBB calculates the effect of the torque 
and the radiation it produces only for radii r > r in = r ms , and does not include the emission 
from smaller radii. The model is thus incomplete. This is one more reason for not including 
a finite torque in fitting spectral data with KERRBB. When 77 = 0, no additional energy 
dissipation is expected for r < r in . In this case, KERRBB includes all the emission and is 
self-consistent. 

We have compared KERRBB with the various other models mentioned above: GRAD, 
DISKPN, EZDISKBB, DISKBB. As §§3.1, 3.2, Figures 10-12 and Tables 1 and 2 show, 
the differences are quite large (except in the case of GRAD), indicating that these models 
make significant errors in computing the spectra of disks. The errors are the result of 
the approximations made by the models, e.g., Newtonian or pseudo- Newtonian dynamics, 
neglect of Doppler and redshift effects, inappropriate torque at the inner edge, etc. When 
the relativistic corrections of Zhang et al. (1997) are applied to DISKBB and EZDISKBB, 
the results are improved somewhat for the low inclination case, but the differences from 
the results of KERRBB are still significant. For the high inclination case the relativistic 
corrections of Zhang et al. (1997) do not improve the nonrelativistic results of DISKBB and 
EZDISKBB at all. 

In §3.2 we showed that, in principle, if we knew M, $ b s and D, then we could determine 
the spin parameter a* of the black hole by fitting spectral data. This was first attempted by 
Zhang et al. (1997) for the black hole X-ray binary J1655 (also, see Sobczak et al. 1999). 
For this approach to work, we need fairly accurate estimates of M, i9 G b s and D (see Table 
3), which are not always available (e.g., the distance to J1550 is very uncertain and so the 
estimates of a* obtained for this source in Fig. 12 are highly suspect). Moreover, even if we 
did have accurate estimates of these parameters, the results would still be uncertain because 
of the questionable blackbody assumption. Non-blackbody effects are usually modeled by 
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means of the spectral hardening factor / co i, and a value / co i = 1.7 is conventionally used 
(Shimura & Takahara 1995). Using this value, and allowing -# obs , M and D to vary over 
their allowed la ranges, we estimate for the source U1543 that a* ~ 0.5 — 0.7 (Table 3). 
If instead we assume / co i ~ 1.5 — 1.6, as in the recent paper of Davis et al. (2004), then 
we find a somewhat more rapidly rotating black hole with a* ~ 0.7 — 0.8. On the other 
hand, if we allow the full range of uncertainty as estimated by Shimura & Takahara (1995), 
viz. / col ~ 1.5 — 1.9, then a* could be as low as ~ 0.3. This illustrates that the biggest 
obstacle to estimating black hole spin via this method is the uncertainty in modeling the 
emitted spectrum from the disk. The crude treatment in terms of a simple hardening factor 
/coi needs to be improved significantly before we can make full use of spectral observations 
of black hole accretion disks. More work along the lines of Davis et al. (2004) is highly 
desirable. 

Other limitations include the fact that the model assumes that the disk is perfectly 
flat with zero vertical thickness, and all the returning radiation is absorbed by the disk and 
reemitted as blackbody radiation. A real disk has a finite thickness which moreover varies 
with radius (i.e., the disk is flared); it probably also has a substantial warp. These com- 
plications will modify the effect of returning radiation and will also change the geometrical 
projection factors needed when computing the observed spectrum. Furthermore, at least a 
part of the returning radiation is scattered by the disk, which adds a high-energy compo- 
nent to the spectrum (Cunningham 1976). Modeling such effects would involve too many 
additional parameters and we have not attempted it. 
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A. Integration of Photon Orbits in a Kerr Spacetime I: The E^ ^ Case 

The general orbit of a photon (indeed for any particle) in a Kerr spacetime is described by 
three constants of motion (Carter 1968): the energy-at-infinity E^, the angular momentum 
about the axis of the black hole L z , and another constant Q (when a = 0, Q + L\ is the 
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square of the total angular momentum). Let us write L z = XE^ and Q = QE^. 3 Then, the 
equations governing the orbital trajectories are separable. Since the system is stationary and 
axisymmetric, only the motions in the r- and ^-directions are required in the calculation of 
the radiation spectrum from the disk. The motion in the r — d plane is governed by (Bardeen 
et al. 1972; Chandrasekhar 1983) 
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dr 
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J tie 



where 
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(Al) 

(A2) 
(A3) 



and r e and d e are the starting values of r and The ± signs in equation (Al) must be 
carefully chosen according to the sign of dr/d'd on the orbit of the photon. 



Define /j, = cos-#, then equation (Al) becomes 

dr , d/j, 



where \i e = cos-# e and 
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For a photon crossing the equatorial plane (which is the case that we are interested in), we 
have Q > 0. Then, both [x\_ and /z 2 are nonnegative. Note, /x+/i 2 = Q/a 2 . 

For a photon emitted by the disk, we have d e = tt/2 and /i e = 0, so fi 2 can never exceed 
H 2 + since otherwise \JW^, and V& become imaginary. Then, the integral over // can be worked 
out with the inverse Jacobian elliptic integral (Byrd & Friedman 1954, eq. 213.00) 
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3 The case of rts is treated separately in Appendix B. 
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where < /i < ji + and 



m 



M /4 + /i 2 _ 



(A8) 



In this paper we use the convention for the modulus in an elliptic integral (e.g., in the 
above inverse Jacobian elliptic integral) which is the same as that of Abramowitz & Stegun 
(1972), and differs from that of Byrd & Friedman (1954) by a square. 

The integral over r can also be worked out with inverse Jacobian elliptic integrals. To 
do so, we need to find out the four roots of R(r) = 0. Since R(r — > ±oo) = oo and 
R(r = 0) = —a 2 Q < 0, R(r) = has at least two real roots. The remaining two roots can be 
either real or complex (in the latter case the two must be complex conjugates of each other). 
We consider the two cases separately: 

Case A. R(r) = has four real roots Let us denote the four roots by r±, r 2 , r 3 , 
and r 4 , with r\ > r 2 > r 3 > r 4 . By equation (A2), we have T\ + r 2 + r 3 + r 4 = and 
^i^2^3^4 = —o?Q < 0; the latter implies that r 4 must be < 0. Physically allowed regions for 
photons are given by R > 0, i.e., r > r\ (region I) and r 3 < r < r 2 (region II). (Although 
R > also for r < r 4 , this region is unphysical since r 4 < 0.) 

In region I (r > ri), the integral over r can be worked out by the following integration 
(Byrd & Friedman 1954, eq. 258.00) 
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(A9) 



where 



m 4 = 



(ri - r 4 )(r 2 - r 3 ) 
(ri - r 3 )(r 2 - r 4 ) 



(A10) 



It is easy to show that 



< 



(r 2 - r 4 )(r - r x ) 
{n - r 4 )(r - r 2 ) 



<1, 



< m 4 < 1 , 



when r\ ^ r 2 . Thus, the inverse Jacobian elliptic function sn 1 in equation (A9) is well 
defined. 
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When T\ = r 2 , the integral in equation (A9) becomes oc sn _1 (l|l) = oo. In this case 
the integral over r needs special treatment. Indeed, when r\ = r 2 , the integral over r can be 
expressed in terms of a logarithm, 



f dr f 
J J 
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where in the last step we have used 2r\ + r 3 + r 4 = 0. 



(All) 



In region II (r 3 < r < r 2 ), the integral over r can be worked out by the following 
integration (Byrd & Friedman 1954, eq. 255.00) 
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It can be shown that 



< 



(ri - r 3 )(r 2 - r) 



< 1 , 



(r 2 - r 3 )(n - r) 

so the inverse Jacobian elliptic function sn -1 in equation (A12) is well defined when r\ ^ r 2 . 
Again, we need special treatment for the case r\ = r 2 : the integral over r is then given 
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Case 5. R(r) = /ias two complex roots and two real roots Let us assume that r\ and 
r 2 are complex, r 3 and r 4 are real and r 3 > r 4 . Then, we must have r\ = r 2 , where * stands 
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for complex conjugate. Since rir 2 r 3 r4 = |ri| 2 r 3 r 4 = —a 2 Q < , r 4 must be < and r 3 must 
be > . [Thus, it is impossible that r 3 = r 4 , since otherwise it gives r 3 = r 4 = 0, which by 
eq. (A2) implies that Q = A — a = and r\ = r 2 = too.] Since R > implies r > r 3 or 
r < r 4 , the physically allowed region for photons is given by r > r 3 . 

Define u = Sft(ri) = 9ft(r 2 ) and t> = ^(ri) = — $5(r 2 ) , where K and S denote the real and 
imaginary parts of a complex number. Then, the integral over r can be worked out with the 
following integration (Byrd & Friedman 1954, eq. 260.00) 

rr dr r r dr 
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It is easy to verify that 



(A - B)r + r 3 B - r 4 A 
(A + B)r - r 3 B - r 4 A 



4AB(r - r 3 )(r - r 4 ) 



< 1 



[(A + B)r - r 3 B - r 4 A] 2 
for r > r 3 . Using the identity (2AB) 2 - [A 2 + B 2 - (r 3 - r 4 ) 2 ] 2 = 4(r 3 - r A ) 2 v 2 > 0, we have 

-2AB <A 2 + B 2 - (r 3 - r 4 ) 2 < 2Afi . 

Then we can show that 

< m 2 < 1 . 

Hence, the inverse Jacobian elliptic function cn -1 in equation (A14) is well defined. 

Therefore, if we can find the four roots of R = 0, we can work out the integral over 
r. The standard procedure for finding the roots of a quartic equation can be found in, e.g., 
Birkhoff & Mac Lane (1965); Zwillinger (2002). 

Each photon emitted by the disk has three possible fates: escaping to infinity, returning 
to the disk, or being captured by the black hole. In the following we discuss the solutions 
for the photon orbits for each of those the possibilities separately. 
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A.l. Solutions for photons escaping to infinity 

We consider an observer at infinity (in practice corresponding to the limit that the 
distance from the black hole to the observer is ^> the radius of the disk) with a polar angle 
$obs, and calculate the orbit of a photon that reaches him/her. As the photon leaves the 
disk, we have d = $ e = ir/2 and /i — pL e — cos-# e = . As the photon reaches the observer, 
we have $ = $obs and \i = /i bs = cos-# bs- 

The position of a photon on the sky as seen by the observer is specified by a pair of 
impact parameters (cx,f3) (see Figure 13 and the second part of Appendix C for details). 
The coordinate a measures the displacement of the photon image perpendicular to the 
projection of the rotation axis of the black hole, and f3 measures the displacement parallel to 
the projection of the axis. The line of sight to the black hole center marks the origin of the 
coordinates, where a — f3 — 0. It can be checked that as the photon reaches the observer, 
on the photon orbit we have dd/dr > (i.e., dfi/dr < 0) if (3 > 0, and dd/dr < (i.e., 
dfi/dr > 0) if /3 < . Therefore, when f3 > 0, the photon must encounter a turning point at 
H — n + : n starts from 0, goes up to then goes down to /z b s (which is < fi + according 
to the earlier argument). When /3 < 0, the photon must not encounter a turning point at 
ji = /!+: \i starts from and monotonically increases to /i bs • 

Based on the above arguments, we can calculate the total integration over /i along the 
path of the photon from the disk to the observer. The results are 



TfS f *L = V° 7«~' v- (A17) 

where J s h stands for the integration along the path of the photon. Note, by definition, r M 
is positive. Substituting equation (A7) into equation (A17), we get 



Tm ^ 04 + pi) 



K(m^) + sign(P) cn 1 



/ /^obs 











(A18) 



where sign(/3) = 1 if f3 > 0, if f3 — 0, and — 1 if f3 < 0; K(m) is the complete elliptic 
integral of the first kind. 

Now let us consider the integration over r. Since the observer is at infinity, the photon 
reaching him/her must have been moving in the allowed region defined by r > r\ when R = 
has four real roots (case A), or the allowed region defined by r > r 3 when R = has two 
complex roots and two real roots (case B). There are then two possibilities for the photon 
during its trip: it has encountered a turning point at r = r + (r + = r\ in case A, r + = r 3 in 
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case B), or it has not encountered any turning point in r. Define 

"°° dr r e dr 



f 



(A19) 



where r e is the radius in the disk where the photon is emitted. Obviously, according to 
equation (A4), a necessary and sufficient condition for the occurrence of a turning point in 
r on the path of the photon is that < r M . Then, the total integration over r along the 
path of the photon from the disk to the observer is 



dr 
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By definition, r^, r e , and r r are all positive. According to equation (A4), we must have 
t~t = Tp. for the orbit of a photon. 

Case A: R = has four real roots When r\ 7^ r 2 , substitute equation (A9) into 
equation (A19) and set r + = r 1 . We obtain 



(A21) 



(A22) 

y/(ri - r 3 )(r 2 - r 4 ) 

Substituting equations (A21) and (A22) into equation (A20) and letting r r = r M , we can 
solve for r e — the radius in the disk where the photon is emitted. The solution is 

r\(r 2 - r 4 ) - r 2 (r 1 - r 4 ) sn 2 ({ 4 |m 4 ) 




r, = 



(r 2 - r 4 ) - (ri - r 4 ) sn 2 (^ 4 |m 4 ) 



where 



6 = - Too) v/(n -r 3 )(r 2 - r 4 ) . 
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Since sn 2 (— ^ 4 |m 4 ) = sn 2 (£ 4 |m 4 ), the solution given by equation (A23) applies whether 
Tfi — is positive or negative, i.e. no matter whether there is a turning point in r or not 
along the path of the photon. Equation (A23) is essentially equivalent to equation (37) of 
Cadez et al. (1998). 

When r\ = r 2 , the integration from r — r± to any r > r\ is infinite (see eq. [All]), so 
= 00 and there cannot be a turning point in r. Then, by equation (All), we have 
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+ In 
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Setting r r = r M , we then obtain the solution for r e 
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where the identity 2ri + r 3 + r 4 = has been used, and 
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The numerical value for r e obtained from equation (A26) should be substituted back in 
equation (A25) (with r r = r M ) to check if r e is a true root, since the square operation in the 
process of solving equation (A25) might produce a false root. 

Case B: R = has two complex roots and two real roots Substituting equation (A14) 
into equation (A19) and setting r + = r 3 , we obtain 
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Substituting equation (A28) and (A29) into equation (A20) and letting r r = r M , we can solve 
for r P . The solution is 



where 



r 4 A - r 3 B - (r 4 A + r 3 B) cn(^ 2 |m 2 ) 
(A-B)- (A + B) cn(f 2 |m 2 ) 



6= (r M -r 00 )v / AB . 



(A30) 
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Since cn(— ^ 2 |m 2 ) = cn(£ 2 |m 2 ), the solution given by equation (A30) applies no matter 
whether r M — r M is positive or negative, i.e. no matter whether there is a turning in r or not 
along the path of the photon. Equation (A30) is essentially equivalent to equation (39) of 
Cadez et al. (1998), although equation (A30) is much simpler. 
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In the above derivation we have implicitly assumed that r + > th, where rn is the radius 
of the horizon of the black hole. If r + < rn, then obviously the orbit of the photon cannot 
have a turning point at r = r + , since if a photon falls into a black hole it cannot get out. 
Therefore, when r + < rn, the solutions given by equations (A23) and (A30) hold only if 
Tn < Too, i-e. when there is no turning in r along the path of the photon. Finally, the 
solution that we have found above for r e represents a physical solution — i.e., it represents a 
photon emitted by the disk — only if it satisfies r in < r e < r out , where r in is the radius of the 
inner boundary of the disk, and r out is the radius of the outer boundary of the disk. 



A photon returning to the disk must satisfy /i = (i.e., d = 7r/2) as it approaches 
the disk at radius r a . Since as it leaves the disk at radius r e it must also have \i = 0, the 
photon must have encountered a turning in d during its trip. Therefore, along the path of 
the photon, /x must first increases from fj, — to fj, — //+, then decreases from fj, — /i + to 
ji = 0. 4 Then, along the path of the photon the total integration over /i is 



The integration over r depends on the allowed region which the photon moves in. In 
case A (R — has four real roots), there are two allowed regions: region I and region II. If 
the photon moves in region I, which is bounded at r = r±, it may encounter a turning point 
at r = 7*1 . If the photon moves in region II, which is bounded at r = r2 and r = r%, it may 
encounter many turning points at r = r2 and r = r^. In case B (R = has two complex 
roots and two real roots), there is only one allowed region that is given by r > r 3 , so the 
photon may encounter a turning point at r = r 3 . In each of these cases the path of the 
photon also depends on the fact whether th appears in the region that the photon moves in, 
and on the direction of the photon as it approaches r = r a in the disk, determined by the 
value of k r = k a (d/dr) a at r = r a , where k a is the four-wavevector of the photon. Let us 
discuss each of these cases separately. 

Case Al: R = has four real roots, region I (r > n) If r± ^ v-i, the fundamental 



4 Of course this assumes that the returning photon crosses the equatorial plane only once. This is reason- 
able since we assume that as the photon reaches the equatorial plane it will be captured by either the disk 
or the gas in the plunging region. 



A. 2. Solutions for photons returning to the disk 




(A32) 
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integral over r is given by equation (A9). Let us define 
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Other relevant integrations are and r e , defined by equations (A21) and (A22) respectively. 

When k r < at r = r a , i.e., the photon approaches r a from the side of r > r a , the 
solution for r e is determined by the equation r M = r r , where the integral over r along the 
path of the photon is (there is no turning point in r) 



r e dr f r e r a \ 

Jr a VR \Jri Jri J 



dr 
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Obviously, if r M > — r a , then r e = oo, which is of course not a solution for a photon 
returning to the disk. So, the solution for r e exists only if < — r a , which is given by 



the solution to r M = r e — r a . The solution is 



where 
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When fc r > at r = r a , i.e. the photon approaches r a from the side of r < r a , the 
integral over r along the path of the photon is 



Jr-e v^R \J ri Jri j ^ a e ' ^ ^ - aJ 



+ r e , 



(r a < < r a + r^) 
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When t m < r a , there is no turning point in r; when r a < r M < r^, there is a turning point 
at r = r±] and when r M > r a + r^, a solution for r e does not exist. Setting r r = r^, we can 
solve for r e . The solution is given by equation (A35) with 6a replaced by 



6a = - r a )^/(n -r 3 )(r 2 - r 4 ) . 



(A38) 



Since sn 2 (— 6a I ^4) — sn 2 (6a|^4), the above solution applies if r M < r a + r^, no matter 
whether there is a turning point in r or not. 
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Of course, when rn > n, there cannot be a turning point in r along the path of the 
photon. Then, when k r > 0, the solution is given by equations (A35) and (A38) if and only 
if < r a 5 

Now let us consider the case r\ = r 2 . In region I the fundamental integral over r is given 
by equation (All). Since the integration from r — r± to any r > r\ is infinite, there cannot 
be a turning point along the path of the photon. Let us define 



1/(7-1 - r 3 )(ri - r 4 ) 



In 



In 



1 + 



2n 



\/{r\ - r 3 )(ri - r 4 

r 2 + r 3 r 4 + 2r x r 



i/(ri -r 3 )(ri - r 4 ) 
\f{r a ~ r 3 )(r a - r 4 ) 



(A39) 



7-1 



(r a - n)^^! -r 3 )(ri - r 4 ) 
1 



(A40) 



In 



VT^ -^3)(^e -r- 4 ) 



i/(ri - r 3 )(ri - r 4 

r 2 + rsr4 _|_ 2nr e 
(r e - ri)i/(ri - r 3 )(ri - r 4 ) 



ri 



We have 



T~a Too 



(A41) 



(A42) 



Then, when k r < at r = r a , the integration over r along the path of the photon is 



R 



= T„ 



(A43) 



The solution for r e is determined by r r = r M . Obviously, if r M > f a — t^, the solution for r e 
does not exist. So, if r M < r a — r^, we have the solution for r e 



r, = 



r 3 r 4 - 


r i 

lan + (rf + r 3 r 4 ) / - r 3 )(ri - r 4 ) 






7„ - 2n/ a/ (ri - r 3 )(n - r 4 ) 


2 


I 



(A44) 



5 A more strict condition should be r M < th, where th = fjj 11 is given by the right hand side of 
eq. (A33) with r a replaced by rn- This condition ensures that r e > t"h- However, since in the final stage of 
our computation we select only the solutions satisfying n n < r < r out , t m < r a and r M < th lead to the same 
final results. (Note, r; n > rn always.) 
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where 



7a = exp 



% - Tp) \f{n - r 3 )(ri - r 4 ) 



(A45) 



When k r > at r = r a , the integration over r along the path of the photon is 

" Ta dr 



Te - T a 



(A46) 



The solution for r e is determined by r r = r M . Since ^= = oo, the solution for r e always 
exists. It is given by equation (A44) with 7 a replaced by 



la = exp 



(r a + r^) ^/(n - r 3 )(ri - r 4 ) 



(A47) 



As in Appendix A.l, the value of r e obtained above for the case of r\ = r 2 should be 
substituted back in the original equation r r = to check if r e is a true solution. 

Case A2: R = has four real roots, region II (r 3 < r < r 2 ) If r*i 7^ r 2 , the fundamental 
integral over r is given by equation (A12). Let us define the following integrals 
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T 32 - T a2 



(A48) 

(A49) 
(A50) 



(A51) 



When rn > t 3 , there cannot be a turning point at r = r 3 . However, there may be a 
turning point at r = r 2 , since r a must be between rn and r 2 . Then, when fc r < at r = r a , 
the integral over r is 



(r a2 < r M < r a2 + r 32 ) 



(A52) 
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When < r a2 , there is no turning point in r; when r a2 < t> < r a2 + r 32 , there is a turning 
point at r = r 2 ; and when r M > r a2 + r 32 , a solution for r e (satisfying r e > r H ) does not exist 
since it would have a turning point at r = r 3 which is impossible. Setting r r = r M , we can 
then solve for r P . The solution is 



r 2 (n 


- r 3 ) 


- n(r 2 - r 3 ) sn 2 (£ 4b |m 4 ) 


(n 


- r 3 ) 


- (r 2 - r 3 ) sn 2 (£ 4fe |m 4 ) 



where 



(A53) 



&& = ^ ( r M - r «2) V^i-rg)^-^) . ( A54) 



Since sn 2 (— £ 4fe |m 4 ) = sn 2 (£ 4fe |m 4 ), the above solution applies if < r a2 + r 32 , no matter 
whether there is a turning point at r = r 2 or not. 

When r H > r 3 and k r > at r = r a , there cannot be a turning point in r along the path 
of the photon. Then, the integral over r is 

r - = f 7r(f -f)7r Tfl - T - (A55) 

Setting r r = r M , we then get the solution for r e . The solution exists only if < r 3a , and is 
given by equation (A53) with £ 4 {, replaced by 

6b = ^(r M + r a2 )-\/(n - r 3 )(r 2 - r 4 ) . (A56) 

When rn < r 3 , the horizon of the black hole does not affect the motion of the photon, 
and the photon moves in a region bounded by r 2 and r 3 . So, the photon can have many 
turning points at r = r 2 and r = r 3 . To properly account for this fact, we define the following 
two numbers 



n = mt\^-j , t-; = r M - 2nr 32 , (A57) 

where 2r 32 is a whole period of the integration over r, and int(a;) means the integer part of 
a real number x. Obviously, we must have < < 2r 32 . 

Then, when rn < r 3 and k r < at r = r a , the integration over r is r r = 2nr 32 + t' t , 
where 



r? _ /r?J ^ = Tfl2 ~ Te2 ' ^ " Ta2) 

(fr r :+fr r :)7n = T a2 + T e2 , (r a2 < rl < T a2 + T 32 ) . (A58) 

^ (r:+/; 2 +/; e )v%=^ 2+27 -3 2 -^ 2 , «>r a2+ r 32 ) 
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When r' < r a2 , there are n turning points at r = r 2 and n turning points at r = r 3 . When 
T a2 < r ^ < r a2 + T 32, there are n + 1 turning points at r = r 2 and n turning points at r = r 3 . 
When > r a2 + r 32 , there are n+1 turning points at r = r 2 and n + 1 turning points at 
r = r 3 . Setting r r = r M , we then obtain the solution for r e . It is given by equation (A53) 
with £ 4 b replaced by 



6 



2V> 



Ta2)y/{ri - r 3 )(r 2 - r 4 ) , 



< r a2 + r 32 ) 

|(r a2 + 2r 32 - r£) v 7 !^"^)!^ 7 ^ , (r^ > r a2 + r 32 ) 



(A59) 



When r H > r 3 and k r > at r = r a , the integration over r is r, 
' {i:-i:)% = r e 2-r a2 , (r><r 3a ) 

< + Irs) VR =T * a + T32 ~ Te2 ' ( T3a < ^ - T3a + T32 ) 

(/£ + Irs + /r?) 7S = T 3« + T 32 + ^2 , « > T 3a + T 32 ) 



2nr 32 + T 1 , where 



T r =\ 



(A60) 



When < T 3a , there are n turning points at r = r 2 and n turning points at r = r 3 . When 
r 3a < < r 3a + r 32 , there are n turning points at r = r 2 and n + 1 turning points at r = r 3 . 
When > r 3a + r 32 , there are n + 1 turning points at r = r 2 and n + 1 turning points at 
r = r 3 . Setting r r = r M , we then obtain the solution for r e . It is given by equation (A53) 
with £ 4 b replaced by 



6 = 



| (t; + r a2 ) ^(n -r 3 )(r 2 -r 4 ) , (t£ < r 3a ) 

| (r 3a + r 32 - t£) v^i - r 3 )(r 2 - r 4 ) , > r 3a ) 



(A61) 



Now let us consider the case of r\ — r 2 . When ri = r 2 , in region II the fundamental 
integral over r is given by equation (A13). Since the integration from r = r\ to r a < r\ = r 2 
is infinite, there cannot be a turning point at r = r 2 . However, there can be a turning point 
at r = r 3 . Let us define 
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(A62) 



+ 



r\ + r 3 r 4 + 2rir„ 



(ri - r a )A/(ri - r 3 )(ri - r 4 ) 
1 

In 



(A63) 



\/( r l - r 3)( r l - r 4) 



V(r e -r 3 )(r e -r 4 ) 
ri - r e 



-38- 



r\ + r 3 r 4 + Znr e 
(ri - r^yjfa - r 3 )(ri - r 4 ) 



We have 



f 



' a dr 



R 



T a -T 3 . 



(A64) 



(A65) 



Then, when k r < at r = r a , the integration over r along the path of the photon is 
(there is no turning point) 

dr 



r e ~ r a 



(A66) 



The solution for r e is determined by r r = r M . Since JJ 1 ^= = oo, the solution for r e always 
exists. It is given equation (A44) with 7 a replaced by 



7o = - exp 



(r a + r M ) a/ - r 3 )(ri - r 4 ) 



(A67) 



When fc r > at r = r„, the integration over r along the path of the photon is 
C % = ^- f ^ (r^<r a - f 3 ) 



TV = 



(A68) 



When r M < f — f 3 , there is no turning point in r; and when > r a — f 3 , there is a 
turning point at r = r 3 . Since ^ = oo, the solution for r e always exists. It is given by 
equation (A44) with 7 a replaced by 



la = 



— exp 

— exp 



% ~ Ty) \/(ri - r 3 )(n - r 4 ) 



(t„ < r a - f 3 ) 



(t" m - f a + 2f 3 ) a/ (ri - r 3 ) (n - r 4 ) , (t„ > f a - f 3 ) 



(A69) 



Of course, when rn > r 3 , there cannot be a turning point. Then, when k r > at r = r a , 
the solution is given by equation (A44) and the first line of equation (A69) (i.e., for the case 
r,, < f a - f 3 only). 

As in Appendix A.l, the value of r e obtained above for the case of r\ = r 2 should be 
substituted back in the original equation r r = to check if r e is a true solution. 

Case 2: R = has two complex roots and two real roots, region r > r 3 The fundamental 
integral over r is given by equation (A14). Let us define 

{A - B)r a + r 3 B - r 4 A 



T a = 



cn 



(A + B)r a - r 3 B - r 4 A 



m 2 



(A70) 
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Other relevant integrations are r M and r e , defined by by equations (A28) and (A29) respec- 
tively. 

When k r < at r = r a , the integral over r along the path of the photon is (there is no 
turning point) 

T ' = £^ = (£-Ovn =T '- T - (A71) 

Obviously, the solution for r e exists only if < — r a . Then, setting r r = r M , we get the 
solution for r e 

_ r 4 A - r 3 B - (r A A + r 3 B) cnfe Q |m 2 ) 
Te ~ (A - B) - (A + B) cn(U\m 2 ) ' 1 j 



where 



6a = (r M + r a )>/IS. (A73) 



When k r > at r = r a , the integral over r along the path of the photon is 



T e , {r„ < T a ) 



(C £) % T ° ■ (A74) 



When r M < r Q , there is no turning point in r; when r a < r M < r^, there is a turning point at 
r = r 3 ; and when r M > r a + , a solution for r e does not exist. Setting r r = r M , we get the 
solution for r e . The solution is given by equation (A72) with £ 2a replaced by 

^ a = (r^-r a )VAB . (A75) 

When r H < r 3 , the above solution applies if r M < r a + r^, no matter whether there is a 
turning point in r or not. When rn > r 3 , there cannot be a turning point at r = r 3 , then 
the above solution applies only if < r a . 



A. 3. Conditions for photons to be captured by the black hole 

To derive the conditions for a photon emitted by the disk to be captured by the black 
hole, we make the following simplified assumption: Any radiation that returns to the equato- 
rial plane inside the inner boundary of the disk is captured by the black hole. This radiation, 
if it is not swallowed by the black hole directly, will be advected or scattered inward by the 
infalling gas, which has a large inward radial velocity (Agol & Krolik 2000). 
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Following Thome (1974), we define a photon capture function C as follow: C — 1 if 
a photon emitted at r e in the disk is eventually captured by the black hole; C — if the 
photon eventually escapes to infinity or returns to the disk. The task in this subsection is 
to find out under what conditions we have C — 1. 

Obviously a sufficient and necessary condition for C — 1 is that the photon neither 
escapes to infinity nor returns to the disk at a radius beyond the inner boundary. In the 
following we translate this condition to mathematical expressions for each case discussed 
in Appendix A. 2. These results are obtained by carefully following the analysis in Ap- 
pendix A. 2. However, for our purpose in this subsection, we follow the orbit of a photon by 
starting from the radius r e where the photon is emitted, rather than the radius at the end 
of the orbit as we did in Appendix A. 2 when discussing the photons returning to the disk. 
Then, k r > at r = r e means that the photon moves to the side r > r e , and k r < at r = r e 
means that the photon moves to the side r < r e . Since the photon must be emitted by the 
disk, we must have r in < r e < r out . 

In the analysis presented below, in each case (Case Al, Case A2, and Case B) the 
symbols have the same meanings as those in the corresponding case in Appendix A. 2. 

Case Al: R = has four real roots, region I (r > r±) When r\ ^ r 2 and r in > r±, let 
us define 



J ri 



dr 
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Then, C — 1 if and only if one of the following conditions is satisfied: 

Condition 1: {r-y < r H < r in ) & (k r < at r = r e ) & (r M > r e — r in ) ; 
Condition 2: (r H < r\ < r in ) & (k r < at r = r e ) & (r e — r in < r M < r e + r ir 



(A76) 



When r\ = r 2 and r in > r\, we define 
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Then, C = 1 if and only if the following condition is satisfied: 
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Condition 3: (r in > r±) & (k r < at r = r e ) & (r M > f in — f e ) . 



Case A2: R = has four real roots, region II (r 3 < r < r 2 ) When r\ ^ r 2 and r in > r 3 , 
let us define 
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m y/{n - r 3 )(r 2 - r 4 ) 
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Then, C = 1 if and only if one of the following conditions is satisfied: 



Condition 4 
Condition 5 
Condition 6 
Condition 7 



(r 3 < r H < r in ) k (k r > at r = r e ) & (r^ > r e2 + r in ) ; 

(r 3 < r H < r in ) & (fc r < at r = r e ) & (r p > r in - r e2 ) ; 

(r H < r 3 < r in ) & (/c r > at r = r e ) & (r e2 + r in < < r e2 + 2r 32 - r in ) ; 

(r H < r 3 < r in ) & (fc r < at r = r e ) k (r in - r e2 < < 2r 32 - r in - r e2 ) . 



When T\ = r 2 and r in > r 3 , we define 
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(ri - rin)^^! - r 3 )(n - r 4 ) 
Then, C = 1 if and only if one of the following conditions is satisfied: 

Condition 8: (r 3 < r H < r- m ) k (k r < at r = r e ) & (r^ > f e — f in ) ; 

Condition 9: (rH < r 3 < r in ) & (k r < at r = r e ) & (f e — f in < < f e + f in — 2f 3 ) . 

Case 2: R = /ias £u;o complex roots and two real roots, region r > r 3 When r in > r 3 , 
let us define 
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cn 



(A - B)r in + r 3 B - r 4 A 
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' ' y/AB [ (A + B)r in - r 3 B - r 4 A 
Then, C = 1 if and only if one of the following conditions is satisfied: 

Condition 10: (r 3 < r H < r in ) k (k r < at r = r e ) & (r M > r e — r in ) ; 
Condition 11: (r H < r 3 < r in ) & (& r < at r = r e ) & (r e — r in < t m < r e + r in ) . 



(A80) 
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B. Integration of Photon Orbits in a Kerr Spacetime II: The rs Case 

A Kerr black hole has an ergosphere inside which the energy of particles and photons can 
be negative (Misner et al. 1973). The radius of the ergosphere is equal to 2M on the equatorial 
plane. When the black hole rotates fast enough, the inner boundary of the disk may enter 
the ergosphere, i.e., r in < 2M. For instance, when r in = r ms , the inner boundary of the disk 
enters the ergosphere of the black hole when a > 0.9428M. When this happens, for some 
photons emitted by the inner part of the disk, — the energy measured at infinity — might 
be close to zero. This needs special treatment since when Eoo = the approach described in 
Appendix A does not apply. Even when E^ is nonzero but the ratio |£ , oo |/£ , / <C 1, where 
Ei is the energy of the photon measured by a local observer at rest with respect to the gas 
at radius tq on the disk, the numerical approach based on the procedure in Appendix A 
cannot give solutions with sufficient precision since then |A|/r and Qjr\ are large numbers. 
The task of this Appendix is to develop special relations for the case of E^ 0. 



Let us define 



^ , kl « i , (bi) 



Em 

and neglect all terms of order e 2 and higher in our calculations. Then we have 

m = R(r)Jh, = -r 2 + 2M(l-^Qr--^-, (B2) 



z 

-2 



Q + L 2 Z V Q' + A'V Q' + A' 2 



where A' = L z /Ei , Q' = Q/E? , which are geometric quantities that do not depend on the 
energy of photons. The motion of the photon in the r — •& plane is then described by 



Jr. 



* if-* (B4) 



The integral over /x can be worked out by the following indefinite integration 



arctan -=== , [i m = \ ~. , x/2 ■ (B5) 



The two roots of R\r) = are r[ — r±(l + 5) and r' 2 = r 2 (l — 5), where 



ri2 = M±,/M 2 --^-, 5= - ~ 2MaX ' £ (B6) 

Q' + X' 2 v/(Q' + A' 2 ) [M 2 (Q> + A' 2 ) - a 2 Q'] 



-43- 



The physically allowed region is given by R > 0, i.e. r_ < r < r+, where r + 
r_ = min (r^, r^). Thus, the integral over r can be worked out by 

<ir f dr 



/R 7 ^ 7 v / ( r +- r ) ( r " r -) 



= arctan 



2r - r\ - 



max(r;,r^), 



(B7) 



Since the physically allowed region is bounded by r + and r_, the photon must return 
to the disk or be captured by the black hole. It cannot escape to infinity. 

First let us consider the case when the photon returns to the disk. Then, we must have 
fi = cos(7r/2) = as the photon approaches the disk at radius r a . Since as it leaves the disk 
at radius r e we also have fi = 0, the photon must have encountered a turning point in $ 
during its trip. Thus, by equation (B5), along the photon orbit the total integral over jj, is 

djj, 



^ = 2 / 
Jo 



TT 



(B8) 



The integral over r along the photon orbit depends on the sign of k r at r = r a — i.e., it 
depends on the direction from which the photon approaches the disk at r a . It can be shown 
that we always have r_ < th for a 2 < M 2 and Q' > 0. Since a photon cannot get out of a 
black hole after it falls in, its orbit cannot have a turning point in r at r = r_. However, it 
can have a turning point at r = r + . Let us define the following integrals 
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(B9) 

(BIO) 
(Bll) 



Since r_ + = r M and r_ < r H < r a < r + , for a photon emitted by the disk and returning 
to the disk, it must have k r < at r = r a (i.e., the photon must approach the disk at r a 
from the side of r > r a ) and it must have encountered a turning point at r = r+. Then, the 
total integral over r is 

dr 



(f+f) 



'R 



e+ 



Setting r r 



TT, we can solve for r e . The solution is given by 

2a\' 



+ r a = r[ + r' 2 = 2M 1 - 



Q + A' 2 



(B12) 



(B13) 
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Of course, a physical solution must satisfy r e + r a > 2r in . According to equation (B13), this 
condition corresponds to 

-2a\'e > - l) (Q' + A' 2 ) . (B14) 

Note that the solution in equation (B13) and the condition in equation (B14) do not depend 
on whether we choose Eiq = E\ a or Eiq = E\ e . 

From the above results it is easy to obtain the condition for the photon to be captured 
by the black hole, which is as follow: either k r < at r = r e , or k r > at r = r e but the 
inequality (B14) is violated. 



C. Calculation of the Constants of Motion 

In the first part of this section we express the constants of motion of a photon (in 
particular, A and Q) in terms of the direction of the velocity of the photon and the radius 
as the photon crosses the equatorial plane. In the second part, we express the constants of 
motion in terms of the impact parameters of the photon at infinity. The former is useful for 
calculating the orbits of the photons that return to the disk or are captured by the black 
hole, while the latter is useful for calculating the orbits of photons that escape to infinity. 
We will also calculate the redshift factor of the photon in each case. 

1. Evaluation of the constants of motion in terms of quantities on the equatorial plane 
We need to relate the local rest frame attached to the disk fluid to the locally nonrotating 
frame. The relation is given by the Lorentz transformation 

e a t = r[e^-^ey , (CI) 

e% = r [-v$, + ey , (C2) 

e a r = e? r)> (C3) 

e a z = 4), (C4) 

where |e", e", e", e"} are the locally nonrotating frame (see, e.g., Bardeen, Press, & Teukolsky 

1972), { e (t)> e ( r )' e (v)' e 0*)} are t ne local rest frame of the disk [e^ = u a is the four-velocity 
of a disk particle], v v is the azimuthal velocity of a disk particle relative to the locally 
nonrotating frame, and T = (l — t>^) is the corresponding Lorentz factor. 

In the coordinates (t,r,ip,z), the locally nonrotating frame are defined by 

X 



d_ 

dip 



(C5) 
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r 2^ V 2 ^ ^ 



A\ 1/2 /9 
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<9r ^ 
1 (d_ 



(C6) 
(C7) 
(C8) 



where A = r 2 - 2Mr + a 2 , A = r 4 + a 2 r(r + 2M), x = (r 2 A/A) 1/2 (the lapse function in the 
equatorial plane), and u; = 2Mar/A (the frame dragging angular velocity in the equatorial 
plane). 

The direction of the velocity of a photon as it crosses the disk plane is specified by 
the normalized four-wavevector of the photon, n a = k a /k^ = k a / 'E h where k a is the four- 
wavevector of the photon, E\ = = —k a e^ = —k a u a is the energy (frequency) of the 
photon measured in the local rest frame. The components of n a in the local rest frame of 
the disk are 



n<«> = 1 , 



= cos 9 , n (r> = sin 9 cos <p , 



(r) 



= sin 9 sin , 



(C9) 



where (9, 0) are spherical coordinates in the local rest frame, with the polar angle 9 measured 
from the normal of the disk, and the azimuthal angle measured from along the disk 
radial direction. 

By equations (C2) and (C9), we have 

k a e« = T [-VyKefo + k a e a W) ] = TE l ( Vip + sin 9 sin 0) . (CIO) 
By equation (C6) and L z = k a (d/d(p) a , we have 
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Equations (CIO) and (Cll) lead to 



r {vy + sin 9 sin 0) E x 



(C12) 



Similarly, equations (CI), (C5), (C9), and E 



-k a {d/dt) a lead to 



Tx(l + sin 9 sin 0)^ + wL 2 



1 L + X \r 2 



1/2 



f w + 



*v + - — 
x V r 



1/2' 



sin 6> sin > , 



(C13) 
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where in the last step equation (C12) has been used. 

Equations (C12) and (C13) can be simplified by introducing L t = r(A/r 2 ) 1 / 2 v ip , = 
Tx + uL\ and VL = u; + x(r 2 /Ay^v^, which are respectively the specific angular momentum, 
the specific energy-at-infinity, and the angular velocity of disk particles. The results are 



L z = 



E n 



1/2 



T sin 9 sin 



E, 



1/2 



TQ sin 9 sin < 



Ex 



(C14) 
(C15) 



Note that, although E\ is always positive, E^ can be negative when the black hole rotates 
so fast that the radius where the photon is emitted is inside the ergosphere of the black hole 
(see Appendix B). Let ip denote the angle between the photon wavevector and the direction 
of efy in the local rest frame. Then cos^ = sin 9 sin 0, and equation (C15) implies that 
Eoo < if and only if ip < ip < ir, where 



= arccos (-^^) ■ 



(C16) 



The solution for ip Q exists for r < 2M , i.e. inside the ergosphere of the black hole. At 
r = 2M (the boundary of the ergosphere), we have ipo = n for a Keplerian disk. 

Equations (C4), (C8), and (C9) lead to 

2 



Q 



d_ 



r 2 (k a e a z ) 2 = r 2 E? cos 2 9 , 



(C17) 



■&=tt/2 



where the first identity comes from the definition of Q (Bardeen et al. 1972). 
From equations (C14), (C15), and (CI 7) we have 

L z V + (A/r 2 ) 1 / 2 Tsm9sm ( f) 



X = 



Q = 



E^ ~ Et + (A/r 2 y/ 2 m sin 9 sin (f) ' 
Q r 2 cos 2 9 

El ~ [Et + (A/r 2 y/ 2 TQ sin 9 sin <f)f 



(C18) 
(C19) 



These two equations relate A and Q to the direction of the photon wavevector (9, 0), the 
radius (r), as well as the quantities specifying the disk motion [L\ E\ Q, and T — of course 
they are not independent) at the point where the orbit of the photon crosses the disk plane. 
For the formulae for calculating L\ E\ Vt ... in the case of a Keplerian disk, see Page & 
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Thorne (1974). The parameters A and Q do not depend on the energy of the photon Ei, 
which is expected under the geometric optics approximation. 

Equation (C18) can be rewritten as 

- Sltf 

nX = 1 ~ E^ + (A/r 2 ) 1/2 ™sm9sm<f) ' (C20) 



which by equation (C15) leads to 



— = . C21 

e, i-nx y ' 



Equation (C21) is useful for calculating the redshift factor of a photon. For example, if a 
photon is emitted from the disk at r = r e and then absorbed by the disk at r = r a , the ratio 
of locally measured energy (frequency) of the photon at r e and r a is simply 

Ea _ e\ - n e L\ i - n a \ 

E le El - Q a Ll 1 - n e A ' ^ ] 

where label "e" means evaluation at r = r e , and label "a" means evaluation at r = r a . 
Equation (C22) implies that if a photon is emitted and returns to the disk at the same 
radius, then the locally measured energy of the photon does not change. 

Given A and Q, we can also calculate the value of cos#. From equations (C19) and 
(C20), we have 

cos0 = ± —y^iM> (C23) 

where the sign of cos 9 should be chosen according to the problem at hand. If the photon is 
leaving the disk, then < 6 < ir/2, and cos^ is positive. If the photon is approaching the 
disk, then n/2 < 9 < n, and cos# is negative. Note, according to equation (C21), 1 — fiA has 
the same sign as E^ since both E t and E^ — Qtf are positive. 6 Since E^ can be negative, 
1 — fiA can also be negative. This is the reason for the absolute value signs in equation (C23). 

2. Evaluation of the constants of motion in terms of photon impact parameters at 
infinity The apparent position of the disk image as seen by an observer is conveniently 
represented by two impact parameters a and (3, measured relative to the direction to the 
center of the black hole. The impact parameters a and (5 are, respectively, the displacement 
of the image perpendicular to the projection of the rotation axis of the black hole on the sky 



6 E jf — fiZ,t = x/r for any particle whose velocity has only an azimuthal component. 
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and the displacement parallel to the projection of the axis. They are related to the conserved 
parameters A and Q by (Cunningham & Bardeen 1973; Cunningham 1975) 

a = -A esc tf obs , P = ± (Q + a 2 cos 2 tf obs - A 2 cot 2 tf obs ) 1/2 , (C24) 

where -# obs is the polar angle of the observer with respect to the rotation axis of the black hole 
(i.e., the inclination angle). From equation (C24) we can solve for A and Q. The solutions 
are 

A = -a sin $ obs , Q = [3 2 + (a 2 - a 2 ) cos 2 $ obs . (C25) 



For ray-tracing, it is more convenient to use polar coordinates in the plane of the disk 
since disk particles are on circular orbits. Figure 13 shows the relation between the impact 
parameters in the observer's sky (the plane S, the image plane) and the polar coordinates 
in the plane of the disk (the plane S'). For a point P on S, labeled by a pair of impact 
parameters (a, (3) , there is a corresponding point P' on S' . P' is obtained from P by drawing 
a straight line from P parallel to the line of sight. The polar coordinates (r', tp') of P' are 
related to the impact parameters (a, f3) by 

a — r' cos tp' , P — r' sin tp' cos -# obs . (C26) 

Therefore, the image of a circular orbit in the disk plane is a circle in S', but an ellipse in S. 
From equations (C25) and (C26) we obtain 

A = — r' cos ip' sin ^ obs , Q = [r 12 — a 2 ) cos 2 i9 obs . (C27) 

These two equations relate the conserved quantities A and Q to the polar coordinates on the 
S' plane which is gridded into many surface elements to make the ray-tracing computation. 7 
The element of solid angle seen by the observer is then 

dad(3 cos $ ohs r'dr' 'dip' 

d"obe - -^2~ = jy 2 , (C28) 

where D is the distance from the observer to the black hole. 



7 The fictitious image plane is gridded as follow: Using the polar coordinates r' and ip' in 5', with the 
origin at the intersection of the line of sight with S' , we divide logr' uniformly from logr' = logn n to 
logr' = log r ou t7 and divide tp uniformly from ip = to ip = 2tt. Because of the light focusing effect of the 
black hole, the image of the inner boundary of the disk on S' has a somewhat larger radius than r ln . So, 
starting the computation from r' = r; n does not miss points in the inner region of the disk. Since r out » r g , 
the image of the outer boundary of the disk on S' has a radius that is almost equal to r out . 
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D. The Radiation Flux of the Disk and the Balance of Energy 

The conservation of angular momentum in a geometrically thin accretion disk rotating 
around a Kerr black hole in the equatorial plane is described by (Novikov & Thorne 1973; 
Page & Thorne 1974) 

l(Mtf- g )=±nrT;\^ H . (Dl) 

Here, as usual, M is the mass accretion rate, IS is the specific angular momentum of disk 
particles, g is the internal torque of the disk which transports angular momentum outward 
in the radial direction, T z is the (p-z component of the stress-energy tensor of the radiation 
produced by the disk, which describes the flow and transport of angular momentum in the 
vertical direction, and H is the half-thickness of the disk. 

Similarly, the conservation of energy in the disk is described by 



d 
dr 



(Mtf -gd)=-A-nrT t \ =H , (D2) 



where is the specific energy of disk particles, Q, is the angular velocity of the disk, and T t z 
is the t-z component of the stress-energy tensor of the radiation, which describes the flow 
and transport of energy in the vertical direction. 

In a local Lorentz frame, the stress-energy tensor of photons can be expressed as an 
integral (Misner et al. 1973; Thorne 1974) 



-tab 



= J n a n b Idtt , (D3) 



where n a = k a /Ei, k a is the four-wavevector of the photon, E\ is the energy of the photon 
measured in the local Lorentz frame (n a does not depend on the energy of the photon), / is 
the intensity of photons which is a function of direction, and dQ is the element of solid angle 
defined by the direction of the three-wavevector of the photon. 

The evaluation of T ab on the surface of the disk can be written in spherical coordinates 
(9, 4>) in the local rest frame of the disk particles, where 9 is measured from the normal of 
the disk surface (see Appendix C). Then we have 

T t z \ z=H = J n t cos 91(9, <f))dQ , T*\ g=H = J n v cos9I(9, 4>)dQ , (D4) 

where dVt = sin 9d9d(f), and we have used n z = = cos 9 (see eq. [C9]). 
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From the definition of n a , we have n t = k t /Ei = —E^/E^ n v = k^/Ei = L z /Ei, where 
-Eqo is the energy- at- infinity of the photon, and L z = XE^ is the angular momentum of the 
photon about the axis of the black hole. Then, by equations (C14) and (C15) we have 

f A\ 1/2 ( A\ 1/2 

n t = -E^-[ — \ rficosV>, n <p = lS+l-£j rcos?/>, (D5) 

where cos^ = sin 9 sin 0, ip is the angle between the wavevector of the photon and the 
direction of e?^ in the local rest frame. 

Substituting equations (D4) and (D5) into equations (Dl) and (D2), we get 

^j-(Ml)-g} = 4vrr(L t F + S) , (D6) 
^- (ME* - gri^ = Anr(E*F + SO.) , (D7) 

where 

F = J cos 9 IdCl (D8) 
is the local energy flux density of the radiation, 

f a\ 1/2 r 

S= — F cos cos 91 dD, (D9) 



is the stress density of the radiation which represents the transport of angular momentum 
in the vertical direction. When S > 0, the angular momentum is transported along the +z 
direction, which means that the radiation exerts a negative torque on the disk. 

Equations (D6) and (D7) have the formal solution (Cunningham 1976; Li 2002a) 

F = F + F s , (D10) 



where 



F = — 

4:71V 



/+ g(4- (ii „4)I(_g (B t_ Mt) -2 i (D11) 

F s = -I(-^(£t _S)Lt)-2 J\e> -SlL^Srdr , (D12) 
the subscript "in" denotes evaluation at the inner boundary of the disk, and 

/ = - Sll))' 2 J\e* - M))^-dr . (D13) 



-51 - 



The quantity g 1VL denotes the value of the internal torque of the disk at the inner boundary. 
The function / has been worked out by Page & Thorne (1974) and is given by their equa- 
tion (15n). Note, F is just the solution for the radiation flux density of the disk when the 
effect of returning radiation is ignored. 

We call the solution given by equation (D10) a "formal solution" since the stress density 
S is a functional of the flux density of the returning radiation in general. 

The flux density F defined by equation (D8) is the "net" flux of energy, which usually 
contains two components: an outgoing component, corresponding to the radiation leaving 
the disk surface; and an ingoing component, corresponding to the radiation approaching the 
disk surface. To decompose the net flux density into an outgoing component and an ingoing 
component, let us write 

/(M) -Un(M), */2<o<« ■ (D14) 

That is, we use 7 out to denote the radiation intensity for the outgoing photons, and I m to 
denote the radiation intensity for the ingoing photons. Then, the outgoing flux density F out 
defined earlier, is given by 

»2w rir/2 



out 



/ cos9I out dtt = / / I out (9,<f>) cos 9 sin 9d9d(f) , (D15) 
Jn+ Jo Jo 



where Q + denotes the whole solid angle above the disk surface, i.e. the solid angle defined 
by0<#<7r/2 and < < 2rr. The ingoing flux density F in is similarly given by 

r 2ir pit 



F m = - cos9I m dtt = - / I in (9,(j)) cos 9 sin 9d9d<j) 

Jtl- Jo Jn/2 

/»27T /»7r/2 

= / / I in (ir- 9, n + cj)) cos 9 sin 9d9d(j) , (D16) 
Jo Jo 



where f2_ denotes the whole solid angle below the disk surface, i.e. the solid angle defined 
by tt/2 < 9 < 7i and < < 2tt. 

In the last line of equation (D16), we have expressed F in as an integral over the reverse 
direction of the wavevector of the incoming photons, i.e. an integral over Q + , which in 
practice is easier to handle. With the above definitions, we have F out > 0, F in > 0, and 
F = F out — F in . Then, by equation (D10), we have 



out 



F + F in + F s . (D17) 



Equation (D17) shows that the returning radiation makes two distinct contributions to the 
outgoing flux density, F in and Fs- F m comes directly from the energy carried by the returning 
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radiation, while Fs represents the work done by the returning radiation on the disk arising 
from the fact that the returning radiation is nonaxisymmetric about the normal of the disk 
surface at the point where the photon crosses the disk (see below). 

The outgoing radiation, which is emitted by the disk locally, is symmetric around the 
normal of the disk surface. However the ingoing radiation, emitted by the disk remotely and 
focused by the gravity of the black hole, must be highly nonsymmetric. Thus we have in the 
frame moving with the disk fluid 

dl jd<f> = , dl-Jdcf) ^ . (D18) 

Because of the presence of cos ip = sin 9 sin <fi in equation (D9) that defines S, only the part 
of 7 that is nonsymmetric around the disk normal contributes to S. By equation (D18), this 
means that only the ingoing component of the radiation contributes to S. Thus we have 

S = (J^J r J 7^(0,0) cos ^ cos 0dft 
1/2 f 

Y 7 in (7r - 0,7r + 0)cos0sin0sin0dfi . (D19) 
Jn+ 



We assume that the radiation emitted by the disk is either isotropic, i.e., 7 out is indepen- 
dent of 9 and 0, or limb-darkened which corresponds to 7 out oc 2 + 3 cos 9 (see Chandrasekhar 
1950; Cunningham 1976). Then, by equation (D15), the outgoing intensity is related to the 
outgoing flux by 

1 I 1 , for isotropic radiation 

7 out = -F out T , T = \ ' 3 , + . • (D20) 

71 [ I + I cos " ' * OY hmb-darkened radiation 

The evaluation of 7 in in a given direction at radius r = r a in the disk can be calculated 
from the outgoing intensity at radius r = r e at the point where the photon comes from. 
Since IsjEf is invariant along the orbit of the photon (Misner et al. 1973), where 7# ; is the 
specific intensity that is related to the intensity 7 by 7 = J °° I El dE[, we have 

7 in (r = r a ) = g 4 I out (r = r e ) , (D21) 

where g = Ei{r = r a )/Ei(r = r e ) is the redshift factor of the returning photon, given in 
equation (C22). 

Since F in and Fs themselves are determined by F out at the radii where the photons 
contributing to F in and F s come from — as equation (D21) shows — equation (D17) must be 
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solved iteratively. First, we ignore the effect of returning radiation and take F out = F . With 
this outgoing flux as an input, we calculate I in by equations (D20) and (D21), then F [n and 
S by equations (D16) and (D19), and F s by equation (D12). Next, we add F in and F s to 
F to obtain a new outgoing flux density F ont , and repeat the calculations to obtain a new 
Fi n and a new Fs- This process is repeated until the solution converges. In practice, five 
iterations are sufficient for convergence. 

The integrals in equations (DI6) and (DI9) can be computed with the ray-tracing 
technique. We divide the solid angle Q + at radius r a into a number of small elements, and 
for each element in a direction defined by (0, <fi) we trace the orbit of the photon — by using 
the formulae provided in Appendixes A. 2 and C — until the orbit crosses the disk at radius 
r e . Then, we calculate 7 out at r e through equation (D20), the redshift factor g through 
equation (C22), and the corresponding I 1D at r a through equation (D21). Thus we obtain 
the integrand in equations (D16) and (D19) for each element of solid angle. Summing the 
contributions gives rise to Fi n and S at r a . We remark that since most of the returning 
radiation comes from the general direction of the black hole, to obtain results with high 
precision it is necessary to divide the solid angle Q + nonuniformly, with many more resolution 
elements in the direction of the black hole. This is most critical when r a 3> M. 

The radiation returning to the disk is, as we have assumed in this paper, absorbed and 
reprocessed by the disk, and then reradiated away. Then, the radiation that permanently 
leaves the disk consists of only two parts: one part escapes to infinity, the other is captured 
by the black hole. The total power of the disk, i.e. the total energy carried by the radiation 
that permanently leaves the disk per unit time as measured by an observer at infinity, is 
(Thorne 1974) 

£tatai = -47r/ T t z \ z=H rdr. (D22) 

J Tin 

Substituting equation (D2) into equation (D22), and using suitable boundary conditions at 
r out (i.e., ps 1 and gVL ps at r out r in ), we obtain £ t otai m terms of M and g ul , which is 
just equation (1) in §2. Equation (1) clearly shows that the power of the disk comes from 
the gravitational binding energy of disk particles and a contribution from a torque at the 
inner boundary of the disk whenever the torque is nonzero. 
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E. The Disk Spectrum as Observed by a Distant Observer 

The specific flux density of the disk radiation as observed by a remote observer, whose 
distance to the black hole is much larger than the size of the disk, is given by 



E n 



J I Eoh dQ obs , (El) 



where I Eobs is the specific intensity of the radiation, E obs is the photon energy, both measured 
by the remote observer, and dQ obs is the element of the solid angle subtended by the image 
of the disk on the observer's sky. 

Using the fact that I El /Ef is invariant along the path of a photon, where E t is the 
energy measured by any local observer on the path (Misner et al. 1973), equation (El) can 
be rewritten as 

F Eobs = J 9 3 I Ecm dQ ohs . (E2) 

E cm is the energy of the photon at its point of emission on the disk as measured by an 
observer located at that point who is corotating with the disk. I Ecin is the specific intensity 
measured by that observer, and 

9 - (E3) 

is the redshift of the photon (see eq. [C21], where = E ohs , E\ = E cm ). 

Suppose the disk radiates like a blackbody. Then the effective temperature of the disk 
measured by a locally corotating observer is simply 



Te ff (r) 



F out {r) 

0"SB 



1/4 



(E4) 



where o"sb is the Stefan-Boltzmann constant, and F out is the outgoing energy flux of the disk 
measured by the locally corotating observer (Appendix D). Suppose the color temperature 
T co i is related to the effective temperature by 

T co \(r) = f co \T cS (r) , (E5) 

where / co i is a constant. Then, the local specific intensity of the radiation emitted by the 
disk is 

2f~ 4 E 3 

j _ J col cm y (Vpi\ 

lE - ~ exp(E em /k B T co i) - 1 ' 



where &b is the Boltzmann constant, T is a function of 9 (the angle between the wavevector 
of the photon emitted by the disk and the normal of the disk surface) and is given by 
equation (D20) when the disk emission is isotropic or limb-darkened. 

Substituting equations (E3), (E5), and (E6) into equation (E2), we have 



The integral in equation (E7) has a form that is suitable for computations with the ray- 
tracing technique (see the second part in App. C). 



We have compared KERRBB with the model GRAD in XSPEC, which considers a 
relativistic disk around a Schwarzschild black hole (Hanawa 1989; Ebisawa et al. 1991, 2003). 
For this comparison, we set a = in KERRBB and ignored the returning radiation and 
limb-darkening since these are the assumptions made by GRAD. (Note that, for the case of 
a Schwarzschild black hole, the effect of returning radiation is not important, see §3.1.) We 
find that, at low inclination angles KERRBB gives results that are consistent with those of 
GRAD. However, at high inclination angles, the results differ. 

We find that the inconsistency between KERRBB and GRAD is caused by the following 
two problems in GRAD: (1) GRAD uses an incorrect formula to calculate the redshift factor 
of a photon approaching the observer (see below). (2) GRAD takes the resolution in the 
azimuthal direction to be independent of the inclination angle. This leads to large errors for 
high inclination angles. Our ray-tracing code uses a resolution proportional to 1/ cos$ b s , 
which gives better numerical accuracy. 

The correct formula for calculating the redshift factor of a photon emitted by a Keplerian 
disk around a Schwarzschild black hole, adapted to the form used in GRAD, is 



where (dip/dt)^ is the angular velocity of the photon evaluated at the point on the disk 
where it is emitted. Equation (Fl) can easily be obtained from equation (C21) by using 
(d(p/dt) p h = AA/r 4 when a = (Bardeen et al. 1972; Chandrasekhar 1983). 

Equation (Fl) differs from equation (B2) in Ebisawa et al. (2003) — which is used in 
GRAD — in the following respects: the power index of (1 — 2M/r) in equation (Fl) is — 1, 




(E7) 



F. Comparison of KERRBB with GRAD 




(Fl) 



-56- 



whereas in equation (B2) in Ebisawa et al. (2003) it is taken to be —0.5. [Note that equa- 
tion (Fl) agrees with equation (A15) in Ebisawa et al. (1991), except that dip p h/ di p h in their 
(A15) should be dip' ph / dtp^ according to their notation.] 

Equation (B2) in Ebisawa et al. (2003) was derived from equations (A6)-(A14) in Ebi- 
sawa et al. (1991). However, equation (A7) in Ebisawa et al. (1991) is wrong; it should 
be 



/ e \ 

Pr 

Pe' 

V *V / LNRO 



^LNRO 



(1 - 2M/ry\dr ph /d9 ph )/(dt ph /d9 ph ) 
;i - 2M/r)- 1/2 r(d6 ph /d6 ph )/(dt ph /d6 ph ) 
1 - 2M/r)- 1 / 2 r(d<p' p Jd9 ph )/(dt ph /d9 ph ) 



(F2) 



where we have used the same notation as in Ebisawa et al. (1991) and we have set G — c — 1. 
Since on a photon orbit in the Schwarzschild spacetime we have ds 2 = = — (1 — 2M/r)dt 2 + 
(l — 2M/r)~ 1 dr 2 +r 2 d6' 2 +r 2 sin 2 6'dip' 2 , it is easy to check that equation (F2) satisfies the null 
condition e 2 — p 2 — p% — p 2 / = (note, = n/2 on the equatorial plane), but equation (A7) 
in Ebisawa et al. (1991) does not. In addition, there are a few other errors and typos in the 
Appendix of Ebisawa et al. (1991): their equation (A8) is incorrect (subsequently corrected 
in Ebisawa et al. 2003), and in the 4x4 matrix in their equation (All), the element on the 
top-right corner should be — Pj, not (5^ (this is obvious since the Lorentz transformation 
matrix should be symmetric). 

If the above correct equations are used, then one can derive a redshift factor that is 
exactly the same as our equation (Fl). 

Unfortunately, equation (B2) of Ebisawa et al. (2003) has been used in the current 
version of GRAD in XSPEC (both XSPEC11 and XSPEC12) at the time the present paper 
was written. In addition, as Ebisawa et al. applied their equation in GRAD, they replaced 
the "— " sign in front of (1 — 2M/r) by a "+" sign. Although the correct formula was 
shown in Ebisawa et al. (1991) — i.e., their equation (A15), except that d(f p h/dt p h should be 
dip'^/dtph — lately it was replaced by an incorrect equation — equation (B2) in Ebisawa et al. 
(2003). 

In Figure 14 we show some examples that compare the results given by KERRBB, those 
given by GRAD, and those given by a modified GRAD where we replaced the following 
statement in GRAD, 



Red = 1.0D0/(1.0D0+SQRT(R0/2.0D0)/(1.0D0-1.0D0/R0)**0.5D0*dphdx/dtdx), 



by 



Red = 1.0D0/(1.0D0-SQRT(R0/2.0D0)/(1.0D0-1.0D0/R0)*dphdx/dtdx) , 
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(i.e., changing the "+" in front of "SQRT" to a and deleting "**0.5"). The parameters 
for the calculated models are: M = 1OM , D = lOkpc, M = 10 19 g sec -1 , / co i = 1, and 
"$obs = 20°, 60°, and 85° in each panel. (As already mentioned, for KERRBB we set a = 0, 
7] = 0, and turned off returning radiation and limb-darkening.) The horizontal axis shows the 
observed photon energy in units of keV. The vertical axis shows the specific photon number 
density calculated by KERRBB (solid lines), GRAD (dashed lines), and the modified GRAD 
(dashed lines), in units of photons keV -1 cm -2 sec -1 . 

Figure 14 shows that, after the incorrect formula in GRAD is replaced with the correct 
one, GRAD agrees quite well with KERRBB. However, for the case of $ bs = 85°, even the 
modified GRAD differs from KERRBB significantly. This is caused by insufficient azimuthal 
resolution in GRAD when the disk is highly inclined. 
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Table 1. Average Ratios of Fitted Parameters for U1543 



Parameter 


KERRBB GRAD DISKPN EZDISKBB DISKBB 


M 
D 


1 1.03 0.61 0.25 0.017 
1 1.01 0.82 0.67 0.31 


Table 2 


. Average Ratios of Fitted Parameters for J 1550 


Parameter 


KERRBB GRAD DISKPN EZDISKBB DISKBB 


M 
D 


1 1.25 1.54 0.76 0.052 
1 1.11 1.47 1.21 0.55 


Table 3. 


Sensitivity of Results on U1543 to Input Parameters 3, 



Adjusted Parameter and Value M (10 18 g/s) a* x 2 P er dof 



tfobs = 22.2° 


2.27 


0.58 


0.87 


tfobs = 19.2° 


2.14 


0.63 


0.88 


M = 1O.4M 


2.00 


0.71 


0.87 


M = 8.4M 


2.43 


0.48 


0.89 


D = 8.0 kpc 


2.65 


0.54 


0.88 


D = 7.0 kpc 


1.81 


0.68 


0.88 


/col = 1.9 


2.66 


0.32 


0.90 


/col = 1-5 


1.76 


0.83 


0.87 


Nominal values b 


2.18 


0.61 


0.88 



a The fits were done with KERRBB using the spectral data obtained 
on MJD 52,467.20 (see Figs. 10 and 12). For each fit, one of the 
following four parameters was varied from its "correct" value by either 
adding or subtracting la: inclination (20.7 ± 1.5 degrees), mass (9.4 
± 1.0 Mq), distance (7.5 ± 0.5 kpc), and spectral hardening factor 
(1.7 ± 0.2). The resulting M, o* and \ 2 for the fit are listed. 

b Values obtained by fixing the parameters at their central values: 
tf obs = 20.7°, M = 9.4M , D = 7.5 kpc, / col = 1.7. 
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0.5 1 1.5 2 

log Radius 

Fig. 1. — The flux density of the radiation returning to the disk for 77 = (upper panel) 
and rj = 00 (lower panel) for a black hole with a = 0.999M. F in is the flux density of the 
returning (ingoing) radiation (eq. [D16]), and F$ is the flux density arising from the angular 
momentum of the returning radiation (eq. [D12]). For comparison, the outgoing flux density 
when the returning radiation is ignored (Fq, eq. [Dll]) is also shown. The disk radius is in 
units of r g = M and the fluxes are in units of 3M e ff/87rr^. Note that F s is negative. 
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Fig. 2. — Photons emitted by the disk consist of three components: some photons return to 
the disk, some go into the black hole, and the rest escape to infinity. The fraction of each 
component in the "total" energy radiated by the disk is shown as a function of the spin of 
the black hole. The definitions of these fractions, which satisfy t rct + ^bh + '•esc = 1, are given 
in eq. (7). Thin lines correspond to the case of a standard Keplerian disk [r] = 0, g m = 0) 
and thick lines correspond to the case of a nonaccreting disk {jj = oo, M = 0, but g in ^ 0). 
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Fig. 3. — The "net" radiation from the disk — i.e., the radiation that permanently leaves the 
disk — consists of two parts: one falls into the black hole, the other escapes to infinity. The 
upper panel shows the fractions of each as a function of the spin of the black hole. The 
definitions of the fractions are given in eq. (9). The lower panel tests the conservation of 
energy, /bh + /esc = 1, in the calculations. It is seen to be satisfied within the errors of the 
code. Thin lines correspond to the case of a standard Keplerian disk (77 = 0) and thick lines 
correspond to the case of a nonaccreting disk (77 = 00). 
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log(l-a/M) 



Fig. 4. — Upper panel: The spin-up function defined in equation (10). The equilibrium spin 
of the black hole [or, the "canonical" state, Thorne (1974)] corresponds to the condition 
d(a/M)/d\nM = 0, which gives a = 0.9983M when the disk radiation is isotropic, and 
a = 0.9986M when the disk radiation is limb-darkened (indicated by the two vertical lines). 
Lower panel: The efficiency of a standard Keplerian disk in converting rest mass into outgoing 
radiation (see the text for definition) as a function of the spin of the black hole. The 
corresponding efficiency when the effect of returning radiation is ignored is shown by the 
dotted line. 
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log E obs 

Fig. 5. — Upper panel: Effect of the spin of the black hole on the observed spectrum of 
the disk. From left to right: a/M = 0,0.5,0.9, and 0.999. Other parameters are: rj = 0, 
$ obs = 30°, M = 1OM , D = lOkpc, M = 10 19 gsec _1 , and / col = 1. Lower panel: Effect 
of the inclination angle of the disk on the observed spectrum. The inclination angles are: 
$obs = 0°,40°,70°, and 85°, as indicated. Other parameters are: 77 = 0, a = 0.9M, M = 
1OM , D = lOkpc, M = 10 19 gsec _1 , and / co i = 1. The energy E ohs is in keV, and the flux 
density N EohB is in units of photons keV -1 cm -2 sec -1 . 



-67- 




log E obs 

Fig. 6. — Effect of the returning radiation on the observed spectrum of an accretion disk. 
The three cases correspond to rj = (upper panel), 77 = 1 (middle panel), and r\ = 00 
(lower panel). The solid line is the spectrum when the returning radiation is included, and 
the dashed line is the spectrum when the returning radiation is ignored. Parameters are: 
a = 0.999M, $ obs = 30°, M = 1OM , D = lOkpc, M cS = 10 19 gsec _1 , and / co i = 1. The 
dotted line in the upper panel (almost coincident with the solid line) is the spectrum when 
the returning radiation is ignored, but M is increased to 1.23 x 10 19 gsec _1 . The dotted line 
in the middle panel is the spectrum when the returning radiation is ignored, and M e g is 
increased to 1.7 x 10 19 gsec _1 . 
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Fig. 7. — Effect of the torque at the inner boundary of the disk. In each panel, the solid 
line is the spectrum of the disk with a torque corresponding to the indicated value of rj (see 
eq. [2]) and the dashed line is the spectrum of a disk with the same value of M e g- but with the 
torque at the inner boundary set to zero. Parameters are: a = (Schwarzschild black hole), 
$ obs = 30°, M = 1OM , D = lOkpc, M eff = 10 19 gsec _1 , and / co i = 1. The dotted lines in the 
last three panels represent the spectra of zero-torque models in which the values of M e g- and 
/coi have been adjusted for the best fit of the corresponding solid lines. The parameters of 
the dotted line models are: M eS = 10 19 gsec _1 and / co i = 1.15 for rj = 0.3; M e n 
and /coi = 1.4 for r] = 1; M eS = 6 x lO^gsec" 1 and / co i = 2.5 for r] = 16.49. 
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Fig. 8. — Similar to Fig. 7 except that a/M = 0.999. The parameters of the dotted line 
models are: M eff = 8.9 x lO^gsec" 1 and / co i = 1.03 for r] = 0.3; M cS = 7.8 x lO^gsec" 1 
and /coi = 1.08 for 77 = 1; M cfr = 5.4 x 10 18 gsec _1 and / co i = 1.25 for i] = 00. 
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log E obs 

Fig. 9. — Effect of limb-darkening on the observed spectrum. The two cases correspond 
to two different disk inclination angles: i? b s = 20° (upper panel) and $ b s = 80° (lower 
panel). The solid line is the spectrum when the disk emission is limb-darkened, and the 
dashed line is the spectrum when the disk emission is isotropic. Other model parameters 
are: rj = 0, a = 0.999M, M = 1OM , D = lOkpc, M = 10 19 gsec -1 , and / co i = 1.5. The 
dotted line in each panel (hardly visible in the upper panel) represents the best fit of the 
solid line with a disk with isotropic emission by adjusting the values of M and / co i. The 
parameters of the dotted lines are: M = 1.17 x 10 19 gsec _1 and / co i = 1.43 in the upper 
panel, M = 0.72 x 10 19 gsec~ 1 and / col = 1.7 in the lower panel. 
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Fig. 10. — Fits of 10 observations of the source U1543 during the high/soft state. The data 
are from Park et al. (2004). Each observation was fitted separately using KERRBB, GRAD, 
DISKPN, EZDISKBB and DISKBB with / col = 1.7. The resulting estimates of the distance 
D and the mass accretion rate M, and the x 2 of the fit, are shown in the three panels. 
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Fig. 11. — Fits of 10 observations of the source J1550 during the high/soft state. The data 
are from Sobczak et al. (2000). Each observation was fitted separately using KERRBB, 
GRAD, DISKPN, EZDISKBB and DISKBB with / col = 1.7. The resulting estimates of the 
distance D and the mass accretion rate M, and the x 2 °f the fit, are shown in the three 
panels. 
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Fig. 12. — Estimates of the black hole spin parameter a/M in U1543, J1550 and J1655. 
The same 10 observations shown in Figs. 10 and 11 were used for the first two sources, 
and two observations from Sobczak et al. (1999) were used for J1655. For each observation, 
the black hole spin was estimated separately using KERRBB, DISKPN, EZDISKBB, and 
DISKBB with / co i = 1.7. No consistent solution was obtained with DISKBB for any of the 
observations of U1543 and J1550. 
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Fig. 13. — Relation between the impact parameters of a photon in the image plane S of an 
observer and polar coordinates in the equatorial plane S' of the black hole. On the image 
plane S, which is perpendicular to the line of sight, each point P is specified by a pair of 
impact parameters {a, (3). In the equatorial plane S', which is perpendicular to the rotation 
axis of the black hole, each point P' is specified by a pair of polar coordinates (r', <//). P 
and P' are connected by a straight line drawn parallel to the direction of the line-of- sight. 
The coordinates (a, (3) of P are related to the coordinates (r',<//) of P' by equation (C26). 
The polar angle of the observer, or equivalently the inclination angle of the disk, is $ bs- 
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Fig. 14. — Comparison of KERRBB with GRAD. The specific photon number density cal- 
culated by KERRBB is shown with solid lines, calculated with GRAD is shown with dashed 
lines, and that calculated with modified GRAD is shown with dashed lines. Three panels 
correspond to three different disk inclination angles as labeled. In the upper and middle 
panels, the dashed lines are almost coincident with the solid lines. 



